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I.  INTRODUCTION 


Objective 

The  objective  of  this  study  is  to  investigate  the  contribution  of  a 
particular  propellant  gas  chemistry  to  the  overall  wear  in  a  gun  tube. 

As  a  continuation  of  two  previous  studiesl«2  conducted  at  Calspan  to 
research  the  pure  carburizing  and  oxidizing  regimes  of  ballistic  chemistry, 
the  present  work  attempts  to  characterize  the  erosion  potential  of  gas 
mixtures  that  represent  the  narrow  band  of  CO/CO2  mole  ratios,  occupied 
by  present  day  large  caliber  gun  propellants.  This  band  of  CO/CO2  ratio 
ranges  from  a  value  of  3,  approximately  the  neutral  point  of  chemical  acti¬ 
vity  for  the  CO/CO2  system,  to  about  19,  the  highest  value  calculated  for 
an  experimental  gun  propellant,  a  low  temperature  RDX  (trimethylene  trini- 
tramine)  compound,  that  was  investigated  during  World  War  II. 3 


Historical  Perspective 

Historically,  the  connection  of  nitramine  propellants  and  high  CO/CO2 
ratios  to  enhanced  gun  tube  wear  probably  dates  from  the  World  War  II 
period.  One  NDRC  study ^  found  that  an  RDX  formulation  was  more  erosive 
than  a  single  base  propellant  (Ml)  with  a  comparably  low  flame  temperature. 
The  experiments  were  conducted  with  both  an  erosion  combustor  fitted  with 
vent  plugs  and  a  caliber  .50  barrel  fixture. 

Other  research5>6  both  during  and  after  the  war,  found  that  pure 
gas  mixtures  of  CO  and  CO2  exhibited  erosivity  at  below  melting  conditions 


1.  C.C.  Morphy  and  E.B.  Fisher,  "The  Role  of  Carburization  in  Gun  Barrel 
Erosion  and  Cracking,"  ARRADCOM  Contractor  Report  ARBRL-CR-00459,  Calspan 
Corporation,  Buffalo,  New  York,  July  1981.  (AD  A102625) 

2.  E.B.  Fisher  and  C.C.  Morphy,  "The  Role  of  Oxygen  in  Gun  Barrel  Erosion 
and  Cracking--A  Shock  Tube  Gun  Investigation,"  ARRADCOM  Contractor  Report 
ARBRL-CR-00427,  Calspan  Corporation,  Buffalo,  New  York,  April  1980.  (AD  A085720) 

3.  F.C.  Kracek,  "Properties  of  Powder  Gas,"  Hypervelocity  Guns  and  The  Con¬ 
trol  of  Gun  Erosion,  Summary  Technical  Report  of  Division  1,  NDRC,  Vol.  1, 
Washington,  D.C.,  1946,  pp.  21-53. 

4.  J.N.  Hobstetter,  "Erosion  of  Gun  Steel  by  Different  Propellants,"  H>~per- 
velocity  Guns  and  The  Control  of  Gun  Erosion,  Summary  Technical  Report 
of  Division  1,  NDRC,  Vol.  1,  Washington,  D.C.,  1946,  pp.  508-328. 

5.  R.C.  Evans,  F.H.  Horn,  Z.M.  Shapiro,  and  R.L.  Wagner,  "The  Chemical 
Erosion  of  Steel  by  Hot  Gases  under  Pressure,"  J.  Physical  and  Colloidal 
Chemistry,  Vol.  51,  1947,  pp.  1404-1429. 

6.  J.C.W.  Frazer,  F.H.  Horn,  and  R.C.  Evans,  "Vent  Plug  Erosion  by  the  CO- 
CO2  Gas  System,"  NDRC  Contractor  Report  A-310,  Johns  Hopkins  University, 
Baltimore,  Md.,  October  1944. 
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in  an  erosion  combustor.  Other  conclusions  from  these  studies;  namely, 
that  erosion  was  independent  of  CO2  concentration  as  long  as  the  CO/CO2 
ratio  was  greater  than  unity  at  equilibrium,  and  that  various  recognized 
carburizing  catalysts  such  as  hydrogen  sulphide,  ammonia  and  free  hydrogen, 
when  added  to  the  test  mixtures,  greatly  increased  vent  plug  erosion,  led 
to  the  further  conclusion  that  carburization  was  playing  an  important 
role  in  thermochemical  erosion  of  gun  barrels.  Several  explanations  of 
this  phemonenon  were  put  forth.  Cracking  of  excess  CO  would  enrich  a  gun 
tube's  surface  with  carbon,  producing  a  drop  in  solidus  temperature  of 
the  steel,  from  1720®K  to  perhaps  1400® K.  The  gun  bore  would,  therefore, 
melt  and  be  lost  at  a  correspondingly  lower  gas  temperature.  Also,  Frazer 
hypothesized  the  formation  and  highly  exothermic  decomposition  of  iron 
penta  carbonyl  early  in  the  ballistic  cycle,  when  tube  surface  temperatures 
were  rising.  The  dissolution  of  Fe  (CO)5  would  augment  the  already 
increasing  surface  heat  flux,  resulting  in  critical  tube  temperatures- being 
prematurely  reached  and  held  over  a  greater  portion  of  the  cycle  time. 

In  recent  years,  the  connection  of  nitramines  and  CO/CO2  ratios  to 
gun  erosion  has  also  been  alluded  to  but  in  a  slightly  different  context. 
Propellant  formulations  containing  the  aforementioned  RDX  and  now,  also, 

HMX  (cyclotetramethylene  tetranitramine)  are  more  often  compared  to  the 
double  base  (nitrocellulose-nitroglycerine)  and  triple  base  (NC-NG-nitro- 
guanidine)  propellants  that  they  have  to  compete  with  for  acceptance  in 
modern,  large  caliber  gun  systems.  The  earlier  comparison  of  low  flame 
temperature  propellants  with  varying  impetus  levels  has  given  way  to  pre¬ 
sent  day  comparisons  of  high  impetus  propellants  with  different  flame 
temperatures.  Of  course,  the  isochoric  flame  temperature  is  still  a  pri¬ 
mary  consideration  when  selecting  a  gun  propellant  since  this  property  is 
directly  linked  to  tube  wear.  Propellants,  such  as  the  nitramines,  that 
possess  both  high  impetus  and  low  flame  temperature  would  appear- to  have 
a  clear  advantage  in  properties.  Unfortunately,  the  recent  research  in 
nitramines  doesn't  reflect  that  clear  advantage. 

Several  laboratory  experiments,  utilizing  a  blow  out  gun^>®»^  have 
indicated  that  nitramines  are  no  more  erosive  than  double  base  propellants 
that  were  comparable  in  either  flame  temperature  or  impetus.  Ballistic 


7.  R.W.  Geene,  J.R.  Ward,  T.L.  Brosseau,  A.  Niiler,  R.  Kirkmire,  and  J.J. 
Rocchio,  "Erosivity  of  a  Nitramine  Propellant,"  Technical  Report  ARBRL-TR- 
02094,  ARRADCOM,  Aberdeen  Proving  Ground,  MD,  August  1978.  (AD  A060590) 

8.  J.R.  Ward  and  R.W.  Geene,  "Erosivity  of  a  Nitramine  Propellant  with  a 
Flame  Temperature  Comparable  to  M30  Propellant,"  Technical  Report  ARBRL- 
MR-02926,  ARRADCOM,  Aberdeen  Proving  Ground,  MD,  June  1979.  (AD  A074346) 

9.  R.  Geene,  B.  Grollman,  A.  Niiler,  A.  Rye,  and  J.R.  Ward,  "Nitramine  Pro¬ 
pellant  Erosivity-III,"  Technical  Report  ARBRL-TR-02278 ,  ARRADCOM,  Aber¬ 
deen  Proving  Ground,  MD,  December  1980.  (AD  A096878) 
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testing  ’  and  tests  performed  with  different  vented  erosion  tc'^ters^' 
have  indicated,  however,  that  nitramine  appears  more  erosive  than  several 
conventional  propellants.  If  nitramines  have  a  distinct  therii;al  advantage 
over  hotter  propellants,  the  conflicting  results  of  contemporar)  studies 
would  indicate  that  under  certain  conditions,  nitramines  are  at  an  eciual].\' 
distinct,  chemical  or  thermochemical  disadvantage.  Recognition  of  this 
chemical  handicap  and  analyzing  its  possible  origin  are  the  purposes  of 
this  work. 

Keeping  in  mind  that  the  chemical  difference  between  nitramine  and 
conventional  propellant  gas  is  more  than  one  of  just  CO/CO i  mole  ratios, 
the  fact  remains  that  it  is  a  major  difference.  Double  base  propellants 
have  low  CO/CQt  ratios  and  nitramines,  like  single  base  propellants  have 
high  CO/CO2  ratios.  This  research  attempts  to  isolate  that  parameter 
from  the  other  thermochemical  phenomena  that  contribute  to  the  complex 
mechanism  of  gun  barrel  erosion,  and  perhaps,  characterize  the  impact  of 
CO/CO2  ratio  on  the  otherwise  desirable  attributes  of  present  and  future 
gun  propellants. 


Method  of  Research 

Research  is  being  performed  by  utilizing  Cal  span's  Shock  Tube  Gun 
(STG)  facility,  which  duplicates  the  thermodynamics  of  gun  chamber  combus¬ 
tion  with  a  tube  chamber  incorporating  a  moving  wall,  i.c.,  the  loading 
face  of  a  gas  driven  piston.  The  gases  compressed  in  the  tube  cliamber 
by  piston  motion  vent  in  a  conventional  manner  by  forcing  a  projectile 
down  the  adjoining  gun  barrel.  During  the  venting  cycle,  the  gases  must 
pass  through  the  nozzled  flow  channel  of  a  gun-steel  sample  positioned  in 
front  of  the  barrel  entrance.  As  a  result,  the  channel  wall  of  the  sample 
experiences  histories  of  pressure,  temperature,  and  forced  convective 
heating,  similar  to  those  of  a  gun  barrel.  Unlike  the  combustion  system 
where  propellant  formulation  dictates  a  set  flame  temperature  and  impetus, 
the  STG  permits  independent  variation  of  any  or  all  conditions  affecting 
ballistic  histories;  and  unlike  smaller  ballistic  compressors  being  utilized 
in  similar  research  work,  the  STG  can  duplicate  the  pressure  cycle  time 
of  most  large  caliber  gun  systems. 


10.  F.A.  Vassallo,  "A  Report  on  The  Erosivity  of  a  Nitramine  Propellant 
When  Fired  in  a  105mm  M68  Tube,"  .ARRADCOM  Contractor  Report  D.AAKIO- 
80-M-3150,  Calspan  Corporation,  Buffalo,  New  York,  February  1081. 

11.  .A.J.  Bracuti,  L.  Bottei,  J.A.  Lannon,  and  L.H.  Caven\’,  "Evaluation  of 
Propellant  Erosivity  with  Vented  Erosion  Apparatus,"  1080  J.AN.NAF 
Propulsion  Meeting,  Vol.  I,  CPIA,  Johns  Hopkins  Universit>-,  Laurel, 

MD,  March  1980. 

12.  F.A.  Vassallo,  "Thermal  and  Erosion  Phenomenology  in  Medium  Caliber 
Anti-Armor  Automatic  Cannons  (MC-AAAC),"  1980  JANNAF  Propulsion  Meeting, 
Vol.  I,  Chemical  Propulsion  Information  .Agency,  Johns  Hopkins  llnivers  it>- , 
Laurel,  MU,  March  1980. 


Several  computer  programs  provide  analytic  support  for  the  STG 
tests.  One  code  models  the  STG  gas  compression  cycle  using  a  van  der  (Vaals 
equation  of  state.  It  computes  gas  pressure  and  temperature,  convective 
heating,  total  heat  input  and  in-wall  temperatures  of  the  test  sample  sur¬ 
face.  This  program  is  used  mainly  in  a  predictive  mode  to  select  gas  compo¬ 
sition  and  peak  pressure  for  an  actual  test.  A  detailed  description  of  the 
code  can  be  found  in  Appendix  A.  Another  code,  described  in  Appendix  B, 
accepts  gas  mixture  and  pressure  data  from  tests  conducted  and  computes 
equilibrium  temperature  and  species  concentration  values  along  isentropes. 
This  information  can  then  be  returned  to  the  preceding  program  to  improve 
its  prediction  of  barrel  heating  and  wear.  STG  testing  and  computer 
analysis,  as  performed  in  this  program,  are  intended  to  expand  present 
knowledge  of  wear  and  erosion  phenomenology. 
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II.  liXPERIMENTAL  PROCEDURES 


Test  Criteria 

As  stated  in  the  introduction,  the  primary  objective  of  this  study 
is  to  quantify  the  erosion  potential  of  propellant  gases  in  terms  of 
equilibrium  CO/CO2  mole  ratio  at  their  respective  flame  temperatures.  To 
accomplish  this  objective,  two  primary  criteria  of  the  test  procedure 
must  be  met . 


Identification  of  Propellant  Gas  CO/CO2  Mole  Ratios 

The  equilibrium  CO/COi  mole  ratios  of  a  representative  group  of  pro¬ 
pellants,  both  conventional  and  experimental,  are  delineated  over  the 
realistic  range  of  ballistic  conditions.  The  mole  ratios  can  be  obtained 
by  Blake  Code  calculation  of  individual  species  concentration  as  a  func¬ 
tion  of  pressure,  temperature,  and  the  equilibrium  constants  of  all  gas 
species  that  make  up  a  given  propellant  gas.  In  this  study,  the  combustion 
gas  concentrations  of  six  propellants  are  used  to  cover  the  realistic 
range  of  equilibrium  CO/COi  mole  ratios  that  exist  at  ballistic  pressures. 
The  propellants  are  illustrated  in  Figure  1.  M8  and  M5  are  conventional 
double  base  propellants.  M50  is  triple-based.  HFP  and  XT  are  nitraraine 
based  and  Ml  is  a  single  base  propellant. 


Modeling  Propellant  Gas  with  Pure  CO  and  CO2 

Within  the  range  of  CO/CO2  ratios  represented  in  Figure  1  by-  the  six 
propellants,  identical  ratios  of  the  two  pure  gases  are  used  in  STG  tests 
to  model  each  propellant.  The  absolute  concentration  of  carbon  monoxide 
and  carbon  dioxide  combined  would  be  about  50  mole  percent  of  actual  pro¬ 
pellant  gas.  The  present  format  of  the  STG  allows  a  maximum  CO-CO2  con¬ 
centration  of  about  45  mole  percent  of  the  mix,  providing  tlie  absolute 
CO2  concentration  does  not  exceed  about  5  mole  percent.  This  limitation 
results  from  disparity  in  mass  between  the  CO  and  C02  molecules.  Increasing 
the  CO2  concentration  lengthens  cycle  time,  which  is  acceptable  if  cx'clc 
time  is  a  test  parameter,  but  can  be  unfavorable  when  attempting  to  com¬ 
pare  tests  of  different  propellant  gas  models.  More  importantly,  increasing 
CO2  concentration  decreases  the  mix  ratio  of  specific  heats  which  results 
in  pressure/temperature  conditions  that  are  not  representative  of  the 
propellant  gas  being  modeled.  That  is,  in  order  to  generate  a  realistic 
"flame  temperature"  in  the  STG,  the  peak  test  pressure  would  have  to  be 
unrealistically  high.  Figure  2  illustrates  the  pure  gas  iscntropcs  generated 
by  the  Isentropic  Equilibrium  Combustion  Code  (described  in  Appendix  B) 
during  adiabatic  compression  in  the  STG.  The  10  percent  CO2  isentrope 
is  shown  as  a  dashed  line  since  it  cannot  be  used  in  the  present  STG  con¬ 
figuration  for  the  reasons  stated  above.  .Also  shown  in  Figure  2  are  the 
six  propellants  whose  gases  are  being  modeled.  Their  nominal  peak  pressure/ 
temperature  conditions  are  indicated. 
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1.  Ec|ui  1  ibrium  pnopcllcinb  gaci  CO/CO—  nabio  versus  pr'essure 


Test  Matrix  Formulation 

Successful  STG  modeling  of  known  propellant  gases  Should  identify  the 
ratio  of  carbon  monoxide  (carburizes  steel  when  cracked  and  absorbed)  to 
carbon  dioxide  (decarburizes  steel  when  absorbed  and  dissociated)  that 
will  yield  the  minimum  amount  of  erosion  and  the  maximum  wear  life  for  a 
gun  barrel  at  characteristic  gas,  and  therefore,  tube  temperatures.  Based 
on  the  above  test  criteria,  the  test  matrix  is  formulated  in  the  following 
ways : 

1.  Each  test  mixture  used  models  the  CO/CO2  ratio  of  one  or  more  known 
propellants  at  their  respective  flame  temperatures,  e.g.,  HFP  (PPL-A-6380) , 
a  nitramine  propellant  (48.6  mole  percent  RDX)  is  represented  by  a  14% 

CO-2. 5%  CO2  mixture  in  STG  tests,  producing  a  characteristic  CO/CO2  ratio 
of  5.6  at  300  MPa,  3300‘'K. 

2.  Each  test  series,  utilizing  the  same  CO/CO2  ratio  includes  runs  that 
represent  the  temperature  conditions  of  all  the  propellants  that  the  parti¬ 
cular  CO/CO2  ratio  models,  e.g.,  a  gas  mixture  of  9%  CO-2. 5%  CO2  with  a 
CO/CO2  ratio  of  3.6  is  tested  at  3000‘’K  (275  MPa)  to  model  M30  (PPL-A-6372) 
and  tested  again  at  3300*K  (275  MPa)  to  model  M5. 

3.  Each  test  series  includes  additional  tests  as  required  to  delineate 
the  threshold  of  erosion  as  a  function  of  gas  temperature  for  a  particular 
gas  mixture. 

4.  In  comparing  the  results  of  tests  incorporating  different  gas  mixtures, 
it  is  beneficial  to  remove  pressure  and  temperature  as  contributing  variables. 
One  method  of  accomplishing  equalization  of  flow  conditions  is  to  maintain 
test  gas  mass  at  a  constant  level,  without  severely  altering  the  ratio  of 
specific  heats. 

The  method  of  maintaining  flow  conditions  that  found  application  in 
the  present  program  is  the  substitution  of  nitrogen  for  carbon  monoxide, 
to  effectively  alter  the  CO/CO2  ratio.  Because  N2  and  CO  have  near  identi¬ 
cal  molecular  weights,  the  total  mass  of  the  gas  mixture  remains  constant. 

The  limiting  assumption  is  one  of  inert  or  near-inert  nitrogen  activity 
under  test  conditions. 

5.  There  are  two  reasons  for  using  two  absolute  CO2  concentrations  in  the 
test  matrix.  First,  the  5  percent  CO2  mixtures  best  model  the  majority 

of  propellant  gases  shown.  However,  the  2.5  percent  CO2  mixture  is  needed 
for  the  highest  CO/CO2  mix  ratio  of  17.2  because  of  the  45.5  mole  percent 
limitation  on  total  CO  and  CO2  in  a  test  mixture.  Second,  in  actual  pro¬ 
pellant  gas,  the  CO2  concentration  has  been  shown  to  exceed  a  limiting 
value  where  an  increase  in  CO2  does  not  produce  a  like  increase  in  erosion. 
However,  a  CO2  concentration  of  5  mole  percent,  or  less,  may  fall  below  the 
limiting  value  for  certain  test  conditions,  where  the  CO2  level  directly 
affects  the  erosion  potential  of  a  gas  mixture.  If  identical  CO/CO2  ratio 
mixtures  with  different  absolute  CO2  concentrations,  i.e.,  2.5  percent 
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and  5  percent,  are  tested  at  similar  flow  conditions,  the  variation  in 
erosion  potential  can  be  quantified  in  terms  of  a  known  change  in  CO2 
levels. 

6.  Surface  conditions  of  the  sample  remain  essentially  frozen  after  a 
test,  since  the  thermochemical  mechanisms  of  interest  proceed  at  a  meaning¬ 
ful  rate,  only  at  elevated  temperatures.  Through  the  use  of  pure  gases 
in  STG  testing,  the  room  temperature  activity  of  propellant  gas  residues 
on  the  steel  surface  is  eliminated. 
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III.  EXPERIMENTAL  MATERIEL 


Shock  Tube  Gun  Facility 

The  ^hock  Tube  Gun,  as  designed  and  developed  by  Calspan  applies  shock 
tube  principles  to  the  study  of  interior  ballistics.  The  facility  con¬ 
sists  of  a  driver  gas  chamber,  a  driven  tube  containing  a  latchable  flying 
piston,  an  instrumented  gas  collection  chamber  and  an  instrumented  gun 
tube  containing  a  projectile. 

One  unique  design  feature  of  the  STG  is  the  driven  piston  which, 
by  its  presence,  affords  physical  isolation  of  the  driver  gas,  normally 
nitrogen,  and  the  chamber  or  test  gas,  which  varies  in  composition  with 
test  objectives.  The  driven  piston,  by  virtue  of  its  mass  in  conjunction 
with  the  projectile's  masi,,  also  controls  the  compression  history  of  the 
test  gas,  enabling  the  facility  to  duplicate  the  interior  ballistic  environ¬ 
ment  of  various  large  caliber  guns,  up  to  and  including  an  8"  howitzer. 

Another  key  design  feature  of  the  STG  is  the  test  section  contained 
within  the  gas  collection  chamber.  It  accepts  specimens  of  various  physical 
configurations  equipped  with  in-wall  thermocouples  and/or  surface  heat  flux 
meters,  and  holds  them  in  place  adjacent  to  the  attached  barrel.  In 
effect,  this  design  offers  a  replaceable,  highly  instrumented  bore  entrance 
to  an  otherwise  conventional  gun  tube. 

With  suitable  variation  in  test  parameters,  made  possible  by  the 
design  of  the  STG,  including  driver  pressure,  piston  mass,  test  gas  composi¬ 
tion,  specimen  configuration,  shot  start  capability  and  projectile  mass, 
one  may  investigate  any  or  all  of  the  internal  factors  affecting. ballistic 
phenomena,  such  as  the  gun  barrel  erosion  and  cracking,  which  the  present 
study  addresses. 


Construction 

Table  1  lists  the  present  structural  dimensions  of  the  STG  as  used 
in  this  study,  and  dictated  by  adiabatic  compression  modeling  and  material 
availability. 

As  shown  in  Figure  3,  the  projectile  launch  components  consist  of  the 
191  mm  driven  tube,  the  instrumented,  high  pressure,  test  gas  collection 
chamber  or  plenum  and  a  30  mm  smooth-bore  barrel.  These  are  supported 
on  a  shock  table  which  is  free  to  move  on  tracks  in  the  direction  of  piston 
motion,  during  the  severe  impulse  loading  caused  by  unbalanced  chamber 
pressure  due  to  test  gas  compression  by  the  decelerating  piston.  This 
floating  mount  system  minimizes  shearing  forces  to  the  supporting  base 
structure,  but  requires  an  adjustable  pneumatic  brake  on  the  driven  tube 
to  absorb  the  axial  loading  on  the  launch  components,  primarily  the  driven 
tube  itself. 
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Table  1.  Shock  Tube  Gun  Characteristics 


Configuration  Data: 


Driven  Tube  I.D. 

0.191  m 

(7.5  in.) 

Driven  Tube  Length 

24.6  m 

(970  in.) 

Piston  Area 

0.0285  m2 

(44.179  in. 2) 

Piston  Mass 

Up  to  91.0 

kg 

(200  Ibm) 

Projectile  Diameter 

30  mm 

(1.181  in.) 

Projectile  Area 

706  mm2 

(1.095  in. 2) 

Projectile  Mass 

Up  to  0.91 

kg 

(2  Ibm) 

Driver  Volume 

0.885  m3 

(54,000  in. 3) 

Chamber  Volume 

2140  mm3 

(130.8  in. 5) 

Pressure  -  at  release  of  projectile 
Barrel  Length 

Variable 
4.57  m 

(180  in. ) 

The  projectile  capture  components  consist  of  a  telescoping  tube 
coupled  to  the  barrel  muzzle,  a  projectile  blast  chamber,  and  a  sand  filled 
tube  to  decelerate  and  catch  the  projectile.  The  telescoping  tube,  the 
purpose  of  which  is  to  permit  independent  motion  of  the  shock  table  and 
blast  chamber,  contains  replaceable  screens  for  measuring  projectile  velo¬ 
city.  The  blast  chamber  reduces  the  noise  and  pressure  levels  as  the 
projectile  exits  the  barrel. 

Figure  4  shows  the  chamber  and  toggle  restraint  system  needed  to 
contain  the  high  chamber  pressures  and  associated  axial  loads.  Chamber 
pressures  are  sensed  using  piezoelectric  transducers.  The  entrance  region 
of  the  launch  tube  can  accommodate  pressure,  heat  flux,  and  erosion  sensing 
devices . 

The  piston,  which  is  used  to  compress  the  test  gas,  is  made  from 
4340  steel  and  weighs  68  kg,  including  the  latching  block  on  its  rear  face 
which  secures  it  at  the  upstream  end  of  the  driven  tube,  prior  to  release 
or  "firing."  Gas  seal  is  maintained  using  "T"  rings  at  the  front  and 
rear  of  the  piston.  Three  brass  wear  rings  or  "bore  riders"  are  used 
to  prevent  steel-to-steel  contact  between  piston  and  tube.  A  buffer  pro¬ 
jection  on  the  face  of  the  piston  and  a  complementary  piston  stop  ring 
at  the  downstream  end  of  the  driven  tube  prevent  direct  impact  of  the 
piston  into  the  test  gas  collection  chamber  in  the  event  that  the  compressed 
test  gas  develops  insufficient  pressure  or  loses  pressure  prematurely  due 
to  chamber  seal  failure. 


Operation 

Operation  of  the  STG  to  collect  test  data  regarding  ballistics, 
heating  and  erosion  follows  a  fixed  experimental  pattern.  Prior  to  a  run 
being  conducted,  the  facility  is  inspected  for  damage  from  the  preceding 
test  and  needed  repairs  are  performed.  Components  such  as  the  driven  tube. 


Figure  4.  Pieton,  chamber,  ond  boggle  reetraint 
componenbe  oF  bbe  Shock  Tube  Gun 


piston,  and  gas  collection  chamber  which  contact  the  test  gas  during  a 
run  are  carefully  cleaned  to  eliminate  contaminants  from  their  surfaces. 
All  wearing  surfaces,  including  piston  "0"  rings  and  bore  riders,  tube 
seals  and  stop  ring  buffers  are  replaced  if  their  wear  limits  are  reached. 
The  piston  is  then  inserted  in  the  upstream  end  of  the  tube  and  latched 
to  the  driver  release  mechanism.  A  new  projectile  is  inserted  in  the 
barrel . 


The  numbered  specimen  is  weighed  prior  to  testing  using  an  analytical 
balance  for  initial  mass,  given  a  final  cleaning  with  freon,  and  then 
installed  in  the  sample  holder  within  the  gas  collection  chamber.  M-11 
mechanical  pressure  gauges  are  also  installed  in  the  chamber.  Thermocouples 
are  inserted  through  the  chamber  wall,  positioned  in  the  sample  wall  and 
then  the  chamber/barrel  assembly  is  lowered  into  position  and  sealed  to 
the  downstream  end  of  the  driven  tube  with  hydraulic  toggles. 

After  installation  of  the  projectile  and  specimen,  the  entire  tube/ 
chamber  cavity  is  evacuated  to  a  pressure  of  2.0  mm  Hg  or  less.  If  vacuum 
is  maintained  for  a  reasonable  time,  indicating  seal  integrity,  the  cavity 
is  purged  with  argon,  re-evacuated,  again  purged  with  argon,  and  evacuated 
for  a  third  time.  The  cavity  is  then  charged  to  the  local  atmospheric 
condition  with  the  required  partial  pressures  of  the  gases  selected  for 
the  te.xt  mixture.  These  partial  pressures  are  dependent  upon  the  mix 
ratio  desired. 

Equations  for  establishing  partial  pressure  settings  are  derived 
from  the  Dalton  model  of  ideal  gas  mixtures,  which  assumes  the  following: 

1.  The  moles  of  mixture,  n,  equals  the  sum  of  the  moles  of  the  component 
gases,  nA  +  n^  +  n(^Q  +  n^Q.,,  where  A,  N,  CO  and  CO2  are  subscripts  referring 
to  argon,  nitrogen,  carbon^monoxide,  and  carbon  dioxide,  respectively. 

2.  Each  component  gas  in  the  mixture  occupies  the  entire  mixture  volume, 

V,  which  in  this  case  is  the  volume  of  the  driven  tube. 

5.  The  temperature,  T,  of  the  components  before  and  after  mixing  remains 
constant . 


4.  The  mixture  pressure,  P,  in  this  case,  1  atmosphere,  is  reasonably 
low,  to  assure  near  ideal  gas  behavior. 


For  the  components; 


PaV  =  n^RT 

PnV  = 


Pco'^  ■  "COPT 

P„„  V  =  n„_  RT 
CO2  CO2 


For  the  mixture:  PV'  =  nRT 
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Since  V/RT  is  a  constant  in  all  the  equations: 


II.  n. 
A  ^ 


*"00 

"co 


PCO 


‘CO 


£  =  £ 
n 


Rewriting; 

P 


[A 

n 


CO 

P 


"co 


_N_ 

P 


PcOt  *^002 


That  is,  for  each  component  of  a  mixture  of  ideal  gases,  the  mole  fraction 
and  the  ratio  of  the  partial  pressure  to  the  total  pressure  are  equal. 

Upon  completion  of  test  gas  charging,  the  mixture  is  given  time  to 
equilibrate  in  the  driven  tube  while  the  required  instrumentation  including 
pressure  transducers  and  thermocouples  are  connected  to  suitable  recording 
devices  and  checked  for  correct  operation.  The  piezoelectric  pressure  output 
is  recorded  by  a  Nicolet  Explorer  III  digital  oscilloscope  and  stored  for 
future  analysis  on  the  scope's  integral  disk  memory.  The  thermocouple 
output  is  recorded  through  an  analog  to  digital  converter  by  a  Rockwell  AIM 
65  microprocessor,  programmed  to  print  out  time  versus  both  temperature 
and  total  heat  input  for  both  data  channels.  Figure  5  is  an  example  of 
plotted  data  taken  from  the  AIM  printer. 

If  the  instrumentation  checks  out  satisfactorily,  the  nitrogen  tank 
farm  valve  is  opened,  the  tube  air  brake  is  charged,  the  driver  is  pres¬ 
surized  to  the  desired  level  for  the  experiment,  recording  devices  are 
activated  and  the  piston  is  released. 

After  exhausting  residual  driver  pressure,  the  air  brake  is  bled, 
and  the  chamber/barrel  assembly  is  decoupled  from  the  tube.  The  specimen 
is  carefully  removed,  inspected,  weighed  and  measured  diametrally  at  four 
axial  locations.  Hard  copy  is  made  of  all  test  data  for  further  reduction 
and  analysis. 


Test  Specimen  Design 

The  primary  objective  of  this  study  is  to  determine  if  certain  propellant 
gas  conditions  enhance  barrel  erosion  and  cracking.  To  correlate  the 
test  data,  the  specimens  used  are  made  nearly  identical  in  shape  and 
composition,  i.e.,  a  4340  steel  cylinder,  38.1  mm  in  length,  31.75  mm  in 
diameter,  and  bored  concentrically  to  12.7  mm,  to  create  a  sonic  flow 
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condition  durins  the  tests.  The  flow  channel  inlet  is  radiused  to  rudinw 
turbulence  and  to  increase  heat  flux  over  a  larger  portion  of  tiic  fioi. 
channel  surface.  The  samples  arc  small  enough  to  fit  in  the  'peciiien  sta-^e 
of  Calspan's  Scanning  lilectron  Microscope  (SEM)  such  that  bore  surface 
examinations  can  be  conducted  without  using  replicas.  ConverselN- ,  t!io 
samples  are  large  enough  to  register  finite  changes  in  weight  and  ''ore 
diameter,  resulting  from  test  conditions.  Mass  changes  arc  measured  in 
tenths  of  a  milligram  on  an  analytical  balance.  Diametral  recess  ion  is 
measured  to  within  10'-*  mm  at  four  specific  axial  locations,  as  si-.o'.-.n  i;; 
Figure  6.  Wso  shown  in  the  figure  are  the  ports  for  in-wall  thernocoup !es 
which  are  used  to  determine  the  integrated  heat  input,  iiimcnsions  or.  ti'.e 
figure  are  given  in  millimeters. 


Figure  6. 


Teal:  epee i men  used  in  i^he  Shock  Tube  Gun 


Heat  Transfer  Instrumentation 

primary  measurement  of  the  study  is  the  amount  of  bore  heating 
associated  with  each  test.  For  this  measurement,  two  in-wall  therihocouples 
are  installed  in  each  sample,  at  distances  approximatelx-  0..'  mm  from  tiie 
hore  surface.  The  method  of  installation  is  shown  in  Figure  ".  I.ich  of 


Figure  7.  lr~wall  thermocouple  iretolled  in 
the  Shook  Tube  Gun  teet  eeotion 


these  thermocouples  independently  may  be  used  to  determine  r.Lt  i',e:ainy 
to  the  bore.  Total  heat  input  is  calculated  from  the  in-wall  t  :ier:;'.ocou;' 1  e '  s 
output  b>’  use  of  methods  de\clopcd  and  reported  i\\'  Ca  1  sji.an .  *  iy'iLfl;,, 
conversion  of  thermocoup le  output  (milln’olts  vs.  tiiaci  to  total  i.ct  ;',oat 
input  per  unit  area  is  accomplished  h>-  using  the  relationshi  iv 


Qft) 

whore  Q(tl 
IT  it) 

k 

C  s 

t 


AT  it)  v-^kc.'t 

is  the  net  bore  heat  input 

is  the  measured  change  in  in-wall  temperature  as 
of  time 

is  tiie  thermal  conduct ivit>- 

is  the  heat  capacity  per  unit  volume 

is  the  time  after  start  of  heating. 


t'unct  1  on 


Data  reduction  procedure  consists  of  calculating  Q(t  i  using  Tcjuation  (1 
at  successive  time  intervals  to  produce  a  curve  of  Qit)  vs.  t  whicl'i  becomes 
asymptotic  to  the  true  heat  input.  To  compute  the  correct  asymptote,  time 
cero  for  the  start  of  heating  must  be  established  accurate!)'  on  the  thermo¬ 
couple  trace.  It  was  concluded  in  a  previous  STG  study-,  tliat  correctly 
shaped  asx'mptotes  were  most  consistent !>•  produced  by  placing  time  tero  at 
the  intersection  of  the  trace  baseline  and  the  average  slope  of  the  initial 
heat  pulse’s  leading  edge. 

Corrected  values  for  Heat  Input  are  taken  from  end  points  on  the 
Q(t)  vs.  t  curves  that  have  been  divided  by  a  correction  factor  to  account 
for  the  heat  flu.x  being  dispersed  over  an  increasing  surface  area  as  it 
passes  through  the  test  specimen's  radial  wall.  This  correction  factor, 

• (r),  based  on  geometric  cons iderations , is  calculated  from  the  relation 

’-,(r)  =  1 

For  the  test  specimen  bore  radius  of  6.53  mm,  ’Ur)  =  .8" 


15.  F.A.  \assaIlo,  "Mathematical  Models  and  Computer  Routines  Used  in 
Evaluation  of  Caseless  .Ammunition  Heat  Transfer,"  Calsj-ian  Rcixort 
G.'I-2948-Z-l ,  Calspan  Corporation,  Buffalo,  New  York,  .lune  19~1  . 


IV.  EXPERIMENTAL  RESULTS 


Gas  Chemistry  Test  Data 

Seven  groups  of  tests  were  conducted  to  quantify  the  thermochemical 
erosion  of  4540  steel  specimens.  The  primary  variable  among  these  groups 
was  CO/CO2  mole  ratio  with  values  of  2.0,  3.6,  5.6,  8.1,  and  17.2  being 
represented.  CO/CO2  ratios  of  2.0  and  8.1  were  each  represented  by  two 
test  groups  containing  2.5  and  5.0  mole  percent  CO2.  Absolute  CO2  concen¬ 
tration,  as  explained  earlier,  was  the  secondary  variable  among  test  groups 
although  groups  with  CO/CO2  mole  ratios  of  3.6,  5.6,  and  17.2  contained 

2.5  percent  CO2  only. 

The  results  of  all  STG  tests  conducted  for  this  program  are  presented 
chronologically  in  Table  2.  Upon  preliminary  inspection,  it  appears  that 
Table  2  contains  many  more  than  seven  test  gas  groups.  However,  the  last 
five  series  containing  an>'where  from  one  to  three  test  shots  are  continua¬ 
tions  of  the  earlier  groups,  and  the  first  test,  i.e.,  95,  which  uses  the 

45.5  percent  CO  "carburizing"  gas  mixture  of  the  previous  program^  was  merely 
intended  as  a  shakedown  run  for  the  STG  facility  and  for  the  new  instrumen¬ 
tation  being  employed  in  this  study. 

The  first  test  group,  consisting  of  tests  96  through  100  and  tests 
130  and  131  employed  a  40.5  percent  CO-5  percent  CO2  mixture  with  an  effec¬ 
tive  CO/CO2  ratio  of  8.1.  This  test  group  models  NT  (PPL-A-6396) ,  a  nitra- 
mine  propellant  containing  44.8  mole  percent  RDX. 

The  second  test  group,  including  tests  101  through  103,  and  127  through 
129,  doubled  the  effective  CO/CO2  ratio  of  the  previous  group  to  17.2  by 
halving  the  CO2  concentration.  The  43  percent  CO-2.5  percent  CO2  mixture 
did  not  model  the  gas  chemistry  of  any  knr  n  propellant  in  the  normal  flame 
temperature  regime,  but  these  tests  did  provide  erosion  data  at  a  finite 
CO/CO2  ratio  that  falls  between  the  "neutral  activity"  point  and  the  full- 
carburizing  level  of  the  previous  program, ^  throughout  the  temperature/ 
pressure  range  of  interest. 

The  third  group,  made  up  of  tests  104  through  106  along  with  tests 
125  and  126,  returned  to  a  CO/CO2  mole  ratio  of  8.1  while  keeping  the  CO2 
concentration  at  2.5  percent.  This  was  accomplished  by  substituting  22.5 
mole  percent  N2  for  a  like  amount  of  CO.  With  respect  to  the  first  test 
group,  the  molecular  concentrations  of  CO  and  CO2  were  both  halved. 

The  fourth  group  of  runs,  tests  107  through  109  and  124  continued  the 
process  of  lowering  the  CO/CO2  ratio  (to  5.6)  by  replacing  more  CO  with 
N2.  This  14  percent  CO- 2. 5  percent  CO2  mixture  in  the  STG  represents  the 
gas  evolved  from  another  nitramine  propellant,  HFP  (PPL-A-6380)  which  con¬ 
tains  48.6  mole  percent  RDX. 
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Table  2.  -las  chemistry  test  data 


I  he  remaining  a'l.-'i'li  of  the  test  gas  mixture  is  argon 
I’a  rent  lies  i  s  indicate  mass  gain  aiul/or  diametral  nroc 


ActiVf*  C0/(;02  I’cak  Teak  Mass**  Diametral  Recession** 

lest  Spec.  Constituent  Ratio  Press.  Tcnip.  l,oss  A  R  C  I) 

No.  No.  l"uj  (’Normalj  (MI’aj  (“K)  (mg)  fmm)  (mm)  (mm)  (mm) 
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The  process  of  approacliing  the  theoretical  neutral  juiint  of  chemical 
activity  was  completed  by  the  fifth  (tests  110-112,  125)  and  sixth  (tests 
115  through  117)  groups  which  further  lowered  the  CO/COj  to  5.0  and  2.(S, 
respective!)’,  while  maintaining  the  CO^  concentration  at  2.5  leicent.  These 
groups  represent  the  equilibrium  CO/COt  ratios  of  most  coinent i ona 1  double 
and  triple-based  propellants. 

The  seventh  test  groups ,  consisting  of  runs  IIS  through  122  duplicates 
its  predecessor's  CO/CO2  ratio  of  2.0  but  at  double  the  COc  content,  i.c., 

5  mole  percent. 

In  addition  to  listing  the  active  constituents  of  each  test  gas  mixture 
Table  2  also  itemizes  peak  STG  chamber  pressure,  the  associated  peak  tempera 
ture  from  the  isentropic  equilibrium  combustion  code  (Appendix  Bi,  mass 
lost  or  gained  by  a  test  specimen  during  its  test  run  and  the  complimentary 
diametral  changes  in  the  specimen's  flow  channel. 

Correlation  of  Erosion  Data  with  Flow  Conditions 

The  flow  conditions  for  an>’  particular  test  in  this  program  are,  to 
a  large  degree,  set  by  the  initial  conditions  of  the  test.  Peak  test  gas 
pressure  is  a  function  of  the  STG  driver  pressure  and  the  mass  of  the  test 
gas  mixture.  The  mass  is  a  constant  throughout  a  group  of  tests  emplO)ing 
the  same  test  gas  mixture  and  a  conscious  effort  was  made  to  minimize  mass 
change  between  test  groups  by  utilizing  the  CO-N2  exchange  method  to  effec¬ 
tively  alter  the  CO/CO2  ratio  while  maintaining  near  constant  mixture  mass. 
Mass  change  teas  unavoidable  when  switching  CO2  concentration  level  hut 
because  CO2  mole  percents  were  always  low;  either  2.5  or  5.i),  the  change 
in  total  mixture  mass  was  limited  to  about  1.1  percent. 

Figure  8  illustrates  the  relationship  between  peak  test  gas  mixture 
pressure  and  specimen  mass  loss.  Also  included  on  the  figure  arc  erosion 
curves  for  oxidizing,  carburizing,  and  inert  (nitrogen)  mixtures  used  in 
the  two  previous  programs. 1^-  While  initially  it  may  appear  that  this  pro¬ 
gram's  mixtures  were  less  active  than  those  of  the  previous  studies,  the 
combination  of  CO  and  CO2  in  one  mixture  produced  equilibrium  ratio.^-  of 
specific  heat  in  the  range  of  1.37  to  1.40,  which  were  measurably  "colder" 
than  values  calculated  for  the  inert  and  45.5  percent  CO  mixtures , typical  1) 
1.41  or  higher.  It  was  therefore  necessary  to  run  the  present  tests  at 
comparative!)’  higher  gas  pressure  in  order  to  generate  the  same  gas  tempera¬ 
tures  of  interest  that  were  sought  in  the  previous  programs.  Figure  S  lias 
been  included  to  point  out  the  apparent  weakness  of  simple  pressure/ erosion 
correlations . 

Figure  9  illustrates  the  relationship  of  specimen  mass  loss  to  the  peak 
test  gas  mixture  temperature.  Because  the  specimen  wall  temiierature ,  at 
which  the  erosion  takes  place,  is  more  closely  linked  to  the  gas  temj'ierature 
than  to  the  gas  pressure.  Figure  9  is  a  better  correlation  of  a  flow  condi¬ 
tion  with  a  surface  reaction  tlian  is  Figure  8.  .As  was  expected,  most  of 
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the  CO/CO2  mixtures  fell  between  the  inert  and  full  carburizing  curves, 
previously  established . 1 ^  Notable  exceptions  are  the  8.1  CO/CO2 ,  5  percent 
CO2  point  ( dai'k  hemicycle)  at  the  40  mg  level  and  the  17.2  CO/CO2,  2.5  per¬ 
cent  point  (hollow  ellipse)  at  the  120  mg  level.  The  former  lies  very  near 
to  the  45.5  percent  CO  curve  and  may  have  fallen  to  the  right  of  it  if  more 
tests  using  this  mixture  could  have  been  performed.  .Alternately,  this 
mixture  also  represents  the  best  compromise  of  potentials  of  any  used  in  the 
program  to  produce  both  low  temperature  carburization  and  high  temperature 
oxidation.  The  latter  point,  mentioned  above,  appears  to  be  substantially 
to  the  right  of  where  it  was  expected  to  fall,  i.e.,  just  to  the  right  of 
the  45.5  percent  CO  curve.  Again,  further  testing  would  probably  have 
clarified  this  discrepancy. 

To  summarize,  the  erosion  potential  of  CO/CO2  mixtures  in  relation  to 
flow  conditions  appears  to  be  buffered  by  the  opposing  chemical  affinities  of 
the  two  molecules.  CO/COo  mixture  ratios  of  from  2.0  to  5.6  produced 
similar  levels  of  erosion,  despite  any  change  in  CO2  concentration.  These 
mixture  ratios  represent  most  of  the  conventional  double  base  and  triple 
base  propellants.  A  CO/CO2  mixture  ratio  of  8.1  appeared  to  be  more  sensi¬ 
tive  to  its  absolute  COt  mole  percent  content.  This  mixture  ratio  is  repre¬ 
sentative  of  the  gas  produced  by  nitramine  base  propellants. 


Correlation  of  Erosion  Data  with  CO/CO2  Ratio 

An  attempt  was  made  to  better  understand  the  flow  conditions/erosion 
data  of  the  previous  figures  from  the  standpoint  of  CO/CO2  ratio  alone. 
Figure  10  represents  the  isotherms  for  the  test  data  in  terms  of  mass  loss 
and  their  respective  CO/CO-,  ratios.  Only  data  from  tests  using  the  primary 
CO-1  concentration,  i.e.,  2.5  mole  percent,  are  shown. 

It  is  apparent  that  overall  erosion  potential  of  the  test  mixtures 
does  not  vary  substantially  with  CO/CO2  ratio,  nor  with  temperature  other 
than  in  a  purely  thermal  way.  However,  the  5.6  CO/CO2  mole  ratio  mixture 
is  the  least  susceptible  to  increasing  temperature,  in  terms  of  its 
erosiveness.  Ironically,  this  mixture  was  used  to  model  HFP,  a  nitramine 
propellant,  whose  normal  flame  temperature  is  too  low  to  take  advantage  of 
this  high  temperature/ low  erosion  characteristic. 


Correlation  of  Erosion  Data  with  COt  Concentration 

Two  of  the  STG  test  gas  mixtures  had  duplicates,  in  terms  of  CO/CO^ 
ratio,  but  with  double  the  amounts  of  CO  and  CO2.  That  is,  a  percent 
CO-2.5  percent  CO2  mixture  and  a  10  percent  CO-5  percent  CO2  mixture  both 
had  a  ratio  of  2.0.  Similarly,  a  20.5  pei'cent  CO-2.5  percent  CO2  mixture 
and  a  40.5  percent  CO-5  percent  CO2  mixture  both  had  a  ratio  of  8.1.  To 
delineate  the  effects  of  temperature  and  CO2  or  CO  concentration  with 
erosivity ,data  from  these  four  test  groups  were  plotted  as  functions  of 
CO/CO^  ratio  and  mass  loss  on  Figure  11. 
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Inunediatc !>•  evident  is  the  inversion  in  erosivity  as  a  thr. vi  ion  oi' 
absolute  CO-CO^  concentration  for  the  two  different  CO  COj  rativi  iiroays. 
Increasing  CO-COj  concentration  for  the  2.0  CO/COj  mi.xture  results  jii  a 
relative  decrease  in  erosion.  Converscl>',  increasing  CO-CO^  coiicent  rat  i on 
for  the  8.1  CO/COj  mi.xture  results  in  a  relative  increase  in  erosion.  This 
would  indicate  that  a  CO/COt  ratio  of  2.0  lies  to  the  oxidizing  side  of 
the  neutral  activity  point  since  a  2  to  1 ,  CO  to  C0-.  reduction  results  in 
increased  erosion,  and  that  the  8.1  CO/COj  ratio  lies  to  the  carburizing 
side  of  the  neutral  point  since  an  8  to  1,  CO  to  CO2  reduction  results  in 
decreased  erosion.  This  finding  is  in  agreement  with  metallurgical  data 
which  places  tlie  neutral  point  of  activity  at  a  C0/C0->  ratio  of  about  .t. 

The  ratio  of  specific  heats  for  the  2:1  mixtures  was  l..T“  and  for  the 
8:1  mixtures,  1.59.  To  a  small  degree,  this  e.xplains  the  grc;iter  oveiaill 
erosion  at  lower  temperature  and  pressure  for  tlie  higher  CO  CO'  ratio 
mixtures  as  shown  on  Figure  11.  Tlieir  higher  total  concentration  of  the 
two  active  gases  would  also  be  a  contributing  factor  to  increased  erosion 
b\’  the  high  CO/CO2  ratio  mixtures.  Additionally,  the  decreasing  disp:ii-it\- 
in  erosit'eness  with  decreasing  temperature  between  mixtures  witii  the  sam.e 
CO/CO2  ratio  but  different  CO^  concent  rat  ions  is  a  limiting  effect  of  Si C 
action  or  cycle  time  (-1  to  0  milliseconds!  on  ciie  erosion  threshold  of  aim- 
particular  gas  mixture  tested. 

What  is  most  noteworthy  is  that  none  of  the  three  factors  mentioned: 
higher  ratio  of  specific  heats,  higher  absolute  CO-COj  content,  and  limited 
STG  cycle  time  in  relation  to  that  of  any  large  caliber  gun,  explain  away 
the  fact  that  at  5450H  the  2.0  CO/CO2  ratio  mixture  had  fallen  close  to 
its  erosion  tlireshold  as  dictated  b)'  the  STG  flow  conditions  init  that  the 
8.1  CO/CO2  ratio  mixture  had  not  reached  a  similar  state  at  5500K  under 
similar  conditions. 


The  importance  of  this  finding  lies  in  its  application  to  guns  of 
varying  calibers  and  action  times,  bespite  the  fact  that  the  tests  conducted 
with  the  STG  in  this  program  were  unable  to  show  the  exact  combination  of 
temperature  and  cycle  time  that  would  pinpoint  the  critical  CO/CO,  propellant 
gas  ratio  in  large  caliber  gun  erosion,  the  data  presented  does  indicate  a 
trend  towards  increasing  chemical  activity  with  increasing  CO  t'02  ratios  at 
comparatively  longer  action  times  and  lower  temperatures.  lliis  finding  is  in 
agreement  with  previous  studies,  where  pure  gas  mixtures  or  propellant  gases 
with  high  CO/CO2  ratios  proved  no  more  erosive  than  low  CO  CO'  ratio  gases 
when  tested  in  combustion  fixtures  with  short  cycle  times  of  1  to  t  millisecond 
but  produced  greater  erosion  when  tested  in  full  scale  gun  fixtures  with 
normal  c\'cle  times  of  over  8  milliseconds. 

Correlation  of  Frosion  Data  with  Heat  I'lux 

Gasdynamic  barrel  erosion  is  the  result  of  lieat  transfer  to  a  melting 
surface,  i.e.,  a  hot-wall  heat  flux.  Ilot-w.all  heat  flux,  in  this  instance. 


is  defined  as  convective  heating  by  the  gas  stream  of  a  surface  at  the 
melting  temperature  of  steel.  The  flow  parameters  of  gas  density,  velocity 
and  temperature  are  included  in  this  heat  flux  calculation.  Therefore, 
this  quantity  may  be  a  better  means  for  correlating  erosion  data  than  either 
pressure  or  temperature  because  the  hot-wall  heat  flux  includes  effects  of 
both.  The  previous  Shock  Tube  Gun  data  correlations  with  temperature  or 
pressure  were  realistic  because  temperature  is  a  function  of  the  pressure 
in  a  ballistic  compressor  during  the  adiabatic  compression  of  similar  gas 
compositions.  However,  for  different  gas  compositions  the  temperature  can 
change  independently  of  the  pressure.  By  using  a  quantity  such  as  hot-wall 
heat  flux,  these  independent  changes  of  temperature  and  pressure  can  be 
taken  into  account. 

To  perform  the  heat  flux  calculation,  flow  conditions  for  each  test 
were  computed,  using  a  combination  of  two  computer  programs  and  the  experi¬ 
mentally  measured  peak  pressures.  The  STG  cycle  model  (.Appendix  .-1)  was  used 
to  create  the  proper  pressure  profile.  .Also  included  in  this  particular 
code  is  a  representation  of  the  flow  through  the  test  sample  and  the  cal¬ 
culation  of  the  convective  heat  flux  using  the  empirical  equation  for 
turbulent  flow  over  a  flat  plate.  The  gas  conditions  of  temperature, 
density  and  velocit)'  are  evaluated  in  this  code.  The  local  hot  wall  heat 
flux  is  computed  and  a  running  summation  of  the  heat  input  to  the  wall  is 
c\alu3ted.  The  computed  total  heat  input  was  compared  with  the  experimental !>■ 
determined  value  for  the  input  gas  mixture.  .A  factor  was  applied  to  the 
heat  flux  calculation  to  bring  the  total  heat  input  into  agreement  with 
the  e.xperimentally  determined  value.  In  this  manner,  the  instantaneous  value 
of  convective  heating,  as  detennined  by  the  flow  conditions  and  exclusive 
of  chemical  heating,  is  believed  to  be  reasonably  correct. 

The  gas  temperature,  as  determined  by  the  STG  model  in  its  current 
state  of  development,  is  only  an  approximate  calculation.  A  more  accurate 
temperature  calculation  is  provided  by  the  equilibrium  combustion  code  that 
is  described  briefly  in  .Appendix  B.  The  code  computes  an  isentropic  com¬ 
pression  of  gases  beginning  with  an  arbitrary  mixture.  In  addition  to 
determining  the  temperature  and  pressure,  the  concentration  of  the  various 
chemical  species  formed  during  the  equilibrium  combustion  process  are  also 
determined . 


.An  approximate  technique  was  used  to  correct  the  heat  flux  calculated 
by  the  STG  model  to  the  more  accurate  temperature  conditions  of  the  equili¬ 
brium  combustion  code.  The  heat  flux  to  a  surface  is  equal  to  the  product 
of  a  coefficient  and  the  temperature  difference  between  the  gas  and  the 
surface , 

(c) 


(4) 


q  =  h  (T^  -  T  J . 

For  the  heat  flux  to  a  flat  plate  in  turbulent  flow,  the  coefficient  is 
functional  !>■  proportional  to  the  gas  density,  velocit\  and  \  isco.s  i  t>' , 
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These  gas  parameters  can  be  expressed  in  terms  of  temperature  as  Uiliov.s: 


P  ~  T"^  through  the  ideal  gas  equation  of  state 
u  '  through  the  energy  equation, 

and  p  '  from  molecular  transport  theory. 


Thus,  the  approximate  dependence  of  h  on  the  gas  tempi. nature  is 


h  ~  [T 
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The  convective  hot  wall  heat  flux  computed  by  the  STd  ci 
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where  Tgs  and  Tj^-g  are  the  respective  peak  values  for  gas  and  sui-face  tempera 
tures.  Similarly,  the  convective  heat  flux  to  a  melt.ng  steel  surface  is 


■‘hw 


T  (T 

ge  ge 


T  ) 
wm 


where  Tgg  is  the  gas  temperature  calculated  by  the  equi  1  ii'vi  ur.  cor.i'U't  ion 
code  and  is  the  solidus  temperature  of  d.'idO  steel,  1'ZiTk,  Hie  corrected 


value  for  the  heat  flux,  ,  is 
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This  represents  an  approximate  value  of  convective  heat  flux  to  a  molting 
surface . 

The  correlation  of  data  from  the  present  program  with  previously  ol'tained 
mass  loss  data^j-  in  terms  of  hot-wall  heat  flux,  is  shown  in  ligure  12. 

The  data  from  the  two  previous  programs  show  the  shift  in  erosion  threshold 
caused  by  insertion  of  one  percent  ox}gcn  in  the  test  gas  in  one  case,  ten 
percent  CO2  in  another  case,  and  45.5  percent  CO  in  a  third  case;  all  witli 
respect  to  the  inert  case  where  pure  nitrogen  was  the  primaip-  test  gas. 

It  should  be  noted  that  in  the  previous  programs  as  well  as  the  current  one. 
the  test  gas  was  composed  of  at  least  54.5  percent  argon  in  order  to  alter 
the  ratio  of  specific  iieats  to  obtain  the  desired  gas  temiieraturc . 

In  reference  1,  it  was  noted  that  a  ratio  of  carbon  monoxide  to  carbon 
dioxide  (CO/CO2)  of  5  was  the  approximate  neutral  point  wlicrc  a  gas  might  be 
expected  to  behave  like  an  inert  gas.  On  Figure  12,  when  the  CO/'CO-  ratio 
was  2,  5.6,  or  5.6,  most  of  the  data  points  fell  along  a  curve  slightl>  to 
the  right  of  the  previously-obtained  inert  curve.  Therefore  the  gas  condi¬ 
tions  represented  by  these  data  points  more  or  less  coincide  uith  tliose  of 
inert  gases.  It  is  noted  that  most  energetic  single-base,  double-base,  and 
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triple-base  gun  propellants  fall  in  this  region.  Two  point- 

during  Runs  118  and  119  with  a  CO/CO-,  ratio  of  2.0,  do  not  fa  1  i  .  i'.  taj  -  oar.t  . 
The  mass  loss  for  Run  118,  as  indicated  on  the  figure,  i  -.va. r 
ttian  e.xpected  and  Run  119  is  somewhat  lower.  Tlie  reasons  f^r  ..,\oar- 

sions  are  currentls'  une.xplained  and  additional  testing  is  j-et.aiiac  •.  oata;;, 
the  true  mass  loss  functions  for  this  test  gas  mi.'cture.  It  siu'uld  .i!.-.'  bi.' 
noted  that  erosion  between  l'  and  li)  milligrams  lies  iti  the  noisi.  .  r, 

the  figure  and  onle  erosion  levels  in  c.xcess  of  tliat  le\el  should  '  oo:.' 1 J 

s i gni f leant , 

The  higher  concentrations  of  CO,  as  indicated  b>'  the  open  and  closed 
hemic\'clic  svanbols  and  the  open  ellipse-like  symbol,  clearl>-  indicate  great;.  ;• 
erosion  for  these  points  tlian  for  the  other  test  gas  ni.xtures  at  thic  sa::;e 
iieating  lei.'el.  Tlie  uppermost  symbols  at  2”()  and  280  d,  man- -see  of  iieat  flux 
show  the  erosion  threshold  is  essentially  m.idway  between  that  of  the  15.8 
percent  CO  and  the  inert  composition.  It  is  noted  tliat  these  points  contain 
40.5  and  15  percent  CO  with  5  and  2.5  percent,  respect i \e 1  ,  of  COp. 

The  data  obtained  during  Run  129  was  presumably  an  outlier  sir.ee  datci 
obtained  at  lesser  conditions  previousl}'  have  e.xhihitcd  si  gni  fi  cant  i>-  iugi'.er 
erosion.  .\t  present  this  low  value  is  unexplained  but  it  can't  be  ignored. 

This  particular  point  is  pivotal  in  that  it  clouds  some  conclusions  tiiat  r.igi'.t 
otlierwise  be  drawn  from  the  data  of  this  program.  The  mass  loss  fro::;  tiiat 
sample  might  be  e.xpected  to  be  three  times  the  level  measured  if  it  were  to 
be  in  agreement  with  the  other  ellipse  shaped  point  at  a  heat  flux  of  2~',' 
J/mm“-scc.  Clearl)',  additional  data  at  this  iiigh  heating  jHnint  siioulj  i^e 
taken  in  order  to  verif>'  or  refute  the  results  of  this  run. 

Other  data  taken  at  a  CO/COj  ratio  of  S.l  but  wi t!i  20.5  percent  CO 
lie  slightl)'  to  the  left  of  the  inert  curve,  showing  les.s  influence  of 
the  higher  C0'C02  ratio  then  the  tests  with  a  higher  concentrat ion  of  CO. 

The  correlation  of  the  test  gas  data  with  convective  hot  wall  heat  flux 
indicate  the  following  trends. 

1.  The  test  witli  gas  mixtures  containing  both  low  levels  of  CO  and  CO^, 
and  a  low  value  of  CO/CO2  ratio  (5.0  to  2),  behaved  quite  similarly  to  pre- 
v'ious  l>'-obtaincd  inert  data.  An  attempt  to  polish  the  surface  of  tiic  test 
specimen  during  this  program  might  well  e.xplain  the  small  shift  from  tb.c 
previousl)'  obtained  inert  gas  curve  (N-.  and  Am 

2.  Higher  levels  of  CO/CO2  ratio  and  ])crcentagc  content  of  C,0  and  COj  wore 
observed  to  shift  the  erosion  threshold  to  lower  values  of  coincctive  heating. 

5.  The  amount  ly  wliich  the  erosion  thrcsliold  was  shifted  a]''i''cars  to  depend 
upon  both  tlic  CO/CO2  ratio  and  on  the  percentage  content  of  CO  and  COj  in 
the  test  gas,  although  none  of  the  data  obtained  on  this  ]>rogram  wliere  COj 
was  incorporated  in  the  gas  mixture  approached  the  level  of  erosion  obtained 
previously  when  45.5  percent  CO  was  tested  alone  with  argon.  The  fact  that 
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increasing  the  amount  of  CO  in  high  CO/COt  ratio  gases  appears  to  decrease 
the  level  of  heating  at  which  the  erosion  threshold  occurs  seems  to  indicate 
a  diffusion  controlled  reaction  where  the  magnitude  of  the  reaction  is 
governed  by  the  amount  of  reactants  diffusing  through  the  surface.  That  is, 
as  the  percentage  of  the  non-equilibrated  reactants  increase,  the  amount  of 
reaction  increases  correspondingly. 

Characterication  of  Test  Specimen  Surfaces 

The  Scanning  Electron  Microscope  was  used  to  help  characterize  the  sur¬ 
face  features  of  each  test  specimen  used  in  this  program.  The  samples  were 
first  sawed  in  half  along  their  flow  axes,  ultrasonical ly  cleaned,  degaussed, 
and  then  photomicrographs  of  a  representative  area  near  t'.e  inlet,  center, 
and  exit  of  each  specimen's  flow  surface  were  made.  Table  5  lists  the  sur¬ 
face  character istics  in  addition  to  pertinent  test  data  for  a  series 
of  specimens  that  represents  the  seven  gas  mixture  groups  tested,  both 
at  near  threshold  and  above  threshold  conditions.  The  photomicrographs  of 
the  fourteen  specimens  tabulated  are  shown  in  Figures  15  through  26. 

As  can  be  seen  in  the  odd  TUimbercd  figures,  near  threshold  flow  condi¬ 
tions  have  minimal  effect  on  the  surface  features  of  the  4540  steel  sjiecimens, 
irrogardless  of  the  gas  mixture  used.  Erosion  is  confined  to  the  leading 
edges  of  larger  machine  marks  that  are  most  exposed  to  convective  flow  and 
least  ca]table  of  dissipating  lieat  other  than  througli  fusion.  The  major 
part  of  eacli  flow  surface  looks  smooth  and  untested  at  all  of  the  locations 
photographed . 

The  oven  numbered  figures  illustrate  the  contrasting  surface  features 
of  seven  tost  specimens  which  experienced  similar  finite  levels  of  erosion 
at  above  threshold  conditions.  The  surface  characteristics  appear  to  be 
dependent  on  flow  channel  location;  i.e.,  localized  flow  factors  sucli  as 
turbulence  level  or  boundary  layer  thickness,  and  local  heat  flux,  that 
effect  tiie  amount  of  melt  and  how  it  is  deposited,  or  if  it  is  deposited 
on  tlie  altered  surface  of  the  steel. 

The  ciianges  in  flow  surface  appearance  seem  also  to  depend  on  the 
cliemistry  of  tlie  gas  that  a  particular  specimen  was  exposed  to  when  tested, 
ihis  finding  is  corroborated  by  comparing  the  photomicrographs  of  the  present 
test  specimens  with  those  from  the  two  previous  reports. 1,2  for  example,  in 
the  "inlet"  photographs  on  Figures  14  and  16,  there  is  evidence  of  the 
clustered  beading  that  was  evident  on  sample  surfaces  eroded  by  oxygen  and 
carbon  dioxide  gas  mixture.s  used  in  the  oxidation  stud>'.-  The  beading  ma>’ 
be  gr.anular  corrosion;  the  preferential  oxidation  of  alpha  iron.  Figures  14 
and  16  represent  the  two  test  groups  with  the  lowest  CO/COs  ratio  and  there¬ 
fore  the  greatest  potential  for  oxidation-type  reactions. 

In  another  comparison  of  present  and  previous  test  specimen  surfaces, 
the  "pebbling"  alluded  to  in  Tabic  .5,  is  prevalent  on  the  "inert”  gas 


specimens  used  in  the  carburization  study.  This  pebbling  is  tii.'igin  to  oo 
intergranular  corrosion,  i.e.,  the  preferential  erosion  of  carbidos  at  tiie 
steel's  grain  boundaries  where  they  arc  concentrated.  "Pebtiling"  is  most 
prominent  when  no  o.xygen,  either  free  or  complexed,  is  present  iti  tlie  test 
atmosphere  and  the  amount  of  pebbling  appears  to  decrease  and  theii  \anish 
as  either  the  CO/COj  ratio  increases  or  the  COj  concentration  increases,  or 
both.  In  short,  "pebbling"  is  not  a  characteristic  feature  produced  b\ 
either  oxidizing  atmospheres  or  those  gases  containing  more  than  about  Ji 
mole  percent  carbon  monoxide,  which  can  be  referred  to  as  carburizing.  in 
the  latter  case,  the  CO  levels  are  probably  high  enough  so  that  the  eutectic 
carbides  within  the  pearlitic  structure  of  the  grain  arc  eroded  in  addit ion 
to  the  intergranular  carbides.  This  may  explain  the  dendritic  formations 
seen  on  sample  surfaces  that  were  exposed  to  CO  concentrations  of  greater 
than  40  mole  percent,  both  in  this  program  and  the  previous  carburization 
study.  Another  feature  shared  by  high  CO  e.xposure  samples  of  these  two 
programs  is  the  "potato  e>e"  or  ivhorl  formations  found  on  the  dounstrea:". 
surfaces  of  these  specimens.  No  explanation  of  this  uni(|ue  surface  charac¬ 
teristic  is  given  at  this  time. 


Tiililc  S.  Characteristics  of  STC  Test  S|us  imeii  I'low  Surfaces 


Flgur«  14.  Surfao*  oKaraot«rl»tlo»  of  •pooimori  no.  77  (tooi.  no.  114)  oftor  59.5  mllllgpom* 
moao  lo.*  In  o  5X  CO  /  ?.  5Z  C02  atmo.phoro  at  322  MPa,  3590* K 


Figur«  15.  Surfae*  oharaot«rl •tio«  of  •pooimon  no.  85  (t««t  .  122)  oftor  noor  nogllglbl* 

mo**  lo**  In  a  1I3Z  CO  /  5Z  C02  abmo*ph*r*  oi.  313  3430*  K 


Figure  16.  Surfao*  oharaoiaris'kio*  of  •p«olm«n  no.  81  (icct  no.  IIB)  of  tor  82.0  milligromi 
mooo  loo*  in  a  10X  CO  /  5Z  C02  otmoophoro  at  325  MPa,  3460* K 


lo**  in  o  9X  CQ  /  2. 5Z  C02  a'tmo»pK«r«  oi  236  HPa,  3310'  K 
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9X  CO  /  2.  5Z  CD2  ai.mo»pK)«r*  at  327  MPa,  3620*  K 


Figure  19.  SurFao*  oh<iroot«rl«tio»  of  •pooimon  no.  71  (twot  no.  108)  of tor  noor  nogllglbl 
moo*  !©••  in  o  14X  CO  /  2. 5X  C02  otmoophoro  ot  203  MPo,  3170" K 


14%  CO  /  2.5%  C02  atmo.ph«r«  at  321  MPa,  3600- K 


Figur«  21.  Curfao*  oharaot«rietio«  of  •p«oiin«n  no.  67  no.  104)  aft«r  n«ar  ncgllgibl 

mas*  loaa  In  a  20. 5X  CO  /  2. 5X  C02  obmosphara  of  245  MPa,  3340 'K 


Figur*  24.  Surfao*  oKaraoiar of  apaoiman  no.  63  (taat  no.  100)  oftar  30.9  mi 1 1 igroma 
maaa  loaa  In  a  40. 5X  CO  /  5X  C02  oimoapKara  at  293  MPa,  3360 *K 


Flgur«  25.  Surfao«  oharciot«rl»tio«  of  •pooimon  no. 66  no.  103)  oftar  n«ar  n«gllglbl 

ma*a  lo*»  In  a  432  CO  /  2.52  C02  otmoaphor*  at  261  MPo,  3382’K 


Surfaoa  oharaoi.aria'tlo*  of  •poolmon  no.  65  (ioot  no.  102)  of  ton  37.7  mllllgr< 
mao*  loo*  in  a  43Z  CO  /  2. 5Z  C02  otmoophor*  at  268  MPo,  3410*K 


V. 


CONCLUSIONS 


The  objective  of  this  program  has  been  to  investigate  the  potential 
of  various  concentrations  and  concentration  ratios  of  carbon  monoxide  and 
carbon  dioxide,  which  make  up  nearly  half  of  all  propellant  gases,  to  erode 
4540  gun  barrel  steel.  As  such,  the  program  deals  exclusively  with  CO-CO^ 
equilibrium,  sometimes  referred  to  as  the  "soot"  equation.  It  does  not 
study  the  "water-gas"  reaction  equilibrium,  nor  the  presence  of  catal\tic 
trace  compounds,  nor  unique  barrel  material  properties;  all  of  which  effect 
the  CO-CO2  activity  in  the  complete  gun  erosion  system.  The  overall  conclu 
sion  from  this  and  previous  studies  would,  therefore,  have  to  be  that  much 
more  research  needs  to  be  carried  out,  particularly  in  regard  to  tlie 
unstudied  mechanisms  cited  above,  in  order  to  complete  the  gun  erosion 
picture . 


In  regard  to  the  present  program  which  dealt  with  erosion  of  4540 
steel  in  varying  CO/CO2  atmospheres,  the  following  conclusions  i>ere  reached 

1.  CO/COt  ratios  of  2.0  to  5.6  e.xhibit  minimal  erosion  above 
the  previously  established  "inert"  level  and  minimal 
dependency  on  CO  or  CO2  concentration  for  their  erosion 
potential  which  is  indicative  of  time  dependent  reactions. 

There  is  some  indication  that,  of  the  three  CO/CO-,  mole 
ratios  used  in  this  program  and  lying  in  the  2.0  to  5.6 
range,  the  5.6  CO/CO2  ratio  is  the  least  susceptible  to 
increasing  erosiveness  with  increasing  temperature.  This 
indication  implies  that  5.6  most  closely  approaches  the 
theoretical  neutral  activity  value  of  all  CO/CO^  ratios 
investigated. 

2.  CO/COt  ratios  of  8.1  and  17.2  shifted  the  erosion  threshold 
to  less  severe  flow  conditions  than  those  of  the  lower 
CO/CO2  ratio  and  previous  "inert”  cases.  Tlicse  iuglier 
CO/COt  ratios  were  also  more  influenced  b>-  the  absolute 
concentrations  of  the  two  gases,  in  terms  of  enhancing 
their  erosion  potential. 
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APPENDIX  A 

SHOCK  TUBE  GUN  COMPUTER  SIMULATION 


Model  Description 

The  mathematical  model  described  in  this  Appendix  simulates  the  opcratioi 
of  the  Shock  Tube  Gun.  The  simulation  provides  a  complete  descrijition  of 
the  entire  cycle  of  the  Shock  Tube  Gun,  beginning  with  release  and  subsequent 
acceleration  of  the  piston  by  the  high  pressure  driver  gas.  As  the  jiiston 
is  accelerated  through  the  driven  tube,  the  simulation  computes  tlic  increase 
in  pressure  and  temperature  of  the  test  gas.  In  addition,  the  simulation 
evaluates  the  total  temperature,  pressure  and  density  in  the  test  gas  collec¬ 
tion  chamber  or  plenum,  and  computes  the  flow  through  the  test  specimen. 

The  heat  flux  to  the  specimen,  the  resulting  temperature  histor>'  at  a  loca¬ 
tion  on  the  surface  of  the  specimen  and  the  temperature  distribution  in  the 
specimen  normal  to  the  surface  are  also  calculated.  Finall}',  the  simulation 
calculates  the  travel  of  a  projectile  through  the  barrel. 

The  objective  of  this  code  is  to  provide  a  means  for  calculating  test 
conditions  for  the  purpose  of  establishing  the  initial  driver  pressure  and 
gas  mixture.  Differences  in  the  ballistic  cycle  of  the  Shock  Tube  Gun  due 
to  gas  composition  are  reflected  through  differences  in  pressure  and  total 
heat  input.  A  primary  use  of  the  code  is  to  help  distinguish  between 
erosion  due  to  melting  and  that  due  to  chemical  effects.  Tins  is  done 
through  computation  of  the  convective  heating  to  the  .'ample  without  chemistr> 
which  provides  a  means  for  comparing  tests  within  a  test  matrix  on  an  equal 
basis  with  respect  to  the  inherent  flow  heating  of  the  test  gas.  Thus, 
excess  material  removal  from  one  gas  mixture  in  comparison  to  another  is 
likely  due  to  chemical  effects.  This,  in  turn,  allows  estimates  to  be  made 
of  the  effective  heat  input  due  to  chemical  effects.  This  can  be  done  b\- 
comparing  the  heat  input  at  the  onset  of  erosion,  or  at  points  of  eiiual 
erosion  between  inert  and  chemically  active  gases. 

The  major  assumption  applied  with  formulating  the  code  was  that  of 
quasi-steady  operation.  That  is,  pressure  waves  and  other  unsteady  aspects 
of  the  event  are  not  calculated.  The  pressure  is  assumed  to  be  constant 
throughout  the  driver  system  and  throughout  the  driven  tube  at  an>-  instant 
of  time  during  the  compression  cycle.  The  individual  gas  constituents  are 
assumed  to  be  mixed  uniformly.  Furthermore,  the  parameters  used  h\-  the  \an 
der  Waals  equation  of  state  for  an  imperfect  gas  and  tlie  temperature  dejien- 
dency  of  specific  heat  are  assumed  to  be  satisfied  b\’  linear  averaging, 
according  to  mole  fraction. 

The  other  limiting  assumption  that  is  currently  employed  in  this  code 
is  that  of  a  frozen  gas  composition;  i.e.,  the  initial  gas  comiios i t i on  is 
assumed  to  be  maintained  throughout  the  ballistic  cycle.  Fhis  assumiition 
influences  the  resulting  temperature  and  pressures  to  some  extent  in  cases 
where  chemical  reactions  and  dissociation  become  important. 
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Piston  Motion 


The  piston  motion  is  evaluated  by  applying  a  force  balance  on  the  piston. 
The  accelerating  force  is  applied  by  the  high-pressure  nitrogen  driver,  500 
to  700  psi,  on  the  upstream  side  of  the  piston.  The  driven  gas  on  the  down¬ 
stream  side  of  the  piston  is  initially  at  atmospheric  pressure.  It  is 
assumed  that  no  gas  leaks  past  the  piston.  This  is  essentially  verified, 
at  least  initially,  by  the  ability  to  evacuate  the  driven  gas  chamber.  The 
nitrogen  gas  is  represented  by  the  ideal  gas  equation  of  state.  The  driver 
conditions  are  calculated  from  conservation  of  energy  principles  which  are 
used  to  continually  evaluate  the  amount  of  energy  that  is  being  transferred 
from  the  gas  to  the  piston.  .As  cited  previously,  this  calculation  is  quasi¬ 
steady  in  that  the  unsteady  expansion  aspects  are  not  considered  and  the 
pressure  is  assumed  to  be  constant  throughout  the  driver  system. 

The  driver  gas  (nitrogen)  properties  are  assumed  to  be  constant  over 
the  range  of  temperature  and  pressure  encountered  during  the  compression 
cycle  and  are  specified  by: 

Equation  of  state  gas  constant,  R  =  55.0  ft/lbf/lbm°R 
Specific  heat  at  constant  volume,  Cv  =  0.177  Btu/lbm-°R 
Specific  heat  at  constant  pressure,  Cp  =  0.248,  and 
Ratio  of  specific  heats,  Y  =  1.4 


The  assumption  of  a  perfect  seal  at  the  piston  infers  the  existence 
of  constant  gas  mass  in  the  driver  system  during  the  cycle  so  that 


I’d  = 


mpRT 


(A-1) 


where  R  and  m  are  constants,  the  driver  volume,  V,  is  expressed  in  terms  of 
piston  travel  and  the  initial  volume  by 


=  'o  "  -S-^p  CA-2) 

where  A.^,  is  the  driver  tube  bore  area. 

Gas  temperature  is  expressed  in  terms  of  internal  energy  where 

E  =  c  m,,T  (A- 5) 

v  D 

and  the  energy  change  resulting  from  work  expended  through  piston  motion  is 
AE  =  p^.A.^Ax^/778.  =  c^m^^AT  (A-4) 
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The  driver  gas  is  initially  at  room  temperature  and  it  is  assu:;,<.d  ti.at  tin 
small  temperature  decrease  during  the  cycle  is  not  influenced  !'■  heat  traa.sf.  r. 

Piston  motion  is  evaluated  by  appl\'ing  a  force  balance  across  the  pisto:. . 
taking  into  account  the  frictional  drag. 

F  =  A^Cpj^  -  P^)  -  P 

where  pp  and  p-p  are  the  respective  pressures  of  the  driver  and  tost  gases, 
and  D  is  the  frictional  drag  of  the  piston  which  is  expressed  in  toi-ms  c-if 
piston  velo''ity  by 

P  = 

Piston  acceleration,  velocity  and  travel  are  then  evaluated. 


u\'  =  a  At 

P  P 

A.\  =  ^^t“  +  V  At 

P  2  o 

Test  Gas  Compression 

The  compression  of  the  test  gas  occurs  as  a  result  of  jiiston  motion. 
Energy  that  is  added  by  virtue  of  the  compression  is  calculated  from  the 
conservation  of  energy  equation.  The  work  done  by  the  piston  on  the  test 
gas  during  this  compression  is  one  term  in  this  conservation  of  energ_\' 
equation.  Other  terms  include  heat  loss  to  the  wall  of  the  tube,  which 
becomes  important  as  the  gas  temperature  rises.  The  other  important  equation 
is  conservation  of  mass.  The  test  chamber  is  not  a  closed  chamber  but  con¬ 
tains  an  e.\haust  port  at  the  downstream  end  where  the  test  sample  is  located. 
Test  gas  is  allowed  to  flow  from  the  test  chamber  through  tlie  test  sample. 
Therefore,  the  mass  in  the  system  is  not  constant  b>’  \irtuc  of  mass  and 
energy  flow  from  the  driven  tube  and  plenum  cliamber  through  the  test  sami'lo. 
Terms  in  the  conservation  of  energx-  and  mass  equations  reflect  this  mass 
and  energy  loss. 

The  equation  of  state  that  applies  to  the  test  gas  is  the  van  dor 
Waals  equation  which  includes  terms  to  e.xprcss  the  nonlinear  relationship 
between  pressure,  densitx'  and  temperature.  The  terms  for  tliis  equation  ai'o 
determined,  as  mentioned  previously,  by  a  linear  averaging  of  the  mole 
fraction  of  the  test  gas  constituents. 

The  test  gas  specific  heat  is  also  assumed  to  bo  for  non-perfect  gas 
and  is  expressed  in  terms  of  a  linear  function  of  tlie  gas  temi'erature .  The 
coefficients  in  this  expression  arc  also  linear  averages  of  the  mole  t' fac¬ 
tion  of  the  test  gas  compositions.  The  test  gas  specific  heat  is  cxin'essed 


(.A-S  I 
l.A-Pl 
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in  terms  of  a  secant  function  in  which  the  product  of  the  specific  heat  and 
the  temperature  yield  the  internal  energy.  This  is  different  from  the  normal 
expression  of  specific  heat,  a  tangent  function,  whereby  the  product  integral 
of  specific  heat  and  temperature  yield  the  internal  energy.  The  technique 
used  here  provides  a  rather  simple  yet  effective  means  for  evaluating  the 
internal  energy  in  a  finite  difference  scheme  with  many  time  steps. 

The  van  der  Waals  equation  of  state  for  the  test  gas; 


P  = 


(A- 10) 


where 


Si 


27 

T 


RT  .  and 

Cl 


V  is  the  specific  volume,  and  Td  and 
pressure  for  the  ith  gas  constituent, 
based  on  the  mole  fraction  of  the  ith 


Pci  are  the  critical  temperature  and 
OL  and  S  are  the  average  quantities 
constituent . 


The  specific  heat  at  constant  volume  is 


S-  '  “S-o  '  S-Ti"  -  -  %)  t 


(A-11) 


where  Cyo  Cyii  are  the  intercept  and  slope  of  the  temperature-dependent 
specific  heat,  Xgi  is  the  mole  fraction  of  the  ith  constituent,  and  ei  is 
the  temperature  exponent  in  the  curve-fitting  relationship  which  includes 
temperature  and  high  density  effects. 

The  test  chamber  mass  balance  is  given  by 

^  =  "Vo  -  (A- 12) 

where  Am  is  the  mass  flow  through  the  test  sample.  Initially  this  flow  is 
assumed  to  be  negligible  and  the  test  chamber  and  barrel  volumes  are  lumped 
together.  When  the  projectile  velocity  exceeds  100  ft/sec,  the  calculation 
of  flow  through  the  test  sample  is  initiated.  This  computation  involves 
determination  whether  sonic  or  subsonic  flow  conditions  exist  within  the 
test  sample.  The  sonic  static  pressure  is  given  by 

p*  =  pi(l  +  ^  (A-15) 

where  Pj  is  the  test  chamber  pressure. 
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If  p*  is  greater  than  the  barrel  pressure,  p2,  Idownstrcam  the  test 

sample)  then  sonic  conditions  exist  and 


At 


(A- 14) 


where  is  the  test  sample  flow  area,  T*  is  the  sonic  static  temj^eratiiro , 
and  At  is  the  computation  time  interval. 


If  subsonic  conditions  exist, 


Am  =  S.u:  A^^((p^ 


""1  1 


(A- 151 


which  is  the  equation  for  flow  through  a  venturi  in  terms  of  the  uj'stream 
and  downstream  line  pressures. 


The  change  in  internal  energy  in  the  chamber  over  the  calculation 
time  interval  is  gi\-en  b)- 


7:’8 


c  AmT.i 
V  1 


(A-lb1 


which  represents  the  compression  work  due  to  the  piston  and  the  loss  in 
enthalpy  due  to  flow  through  the  test  sample.  Heat  transfer  losses  in  the 
driven  tube  and  chamber  are  not  included  in  this  anal}sis  at  present.  The 
gas  chamber  temperature  is  defined  by 

”^1  ^  (A-n 


where  Ei  is  the  current  value  of  internal  energy. 


Test  Sample  Flow 

The  test  sample  contains  a  straight  cylindrical  cliannel  with  a  raJiuseJ 
entrance.  The  flow  through  the  channel  is  computed  by  cither  sonic  or  sub¬ 
sonic  conditions  depending  on  the  pressure  at  the  inlet  and  outlet  to  the 
sample.  Code  calculations  provide  the  static  flow  conditions  of  pressure, 
temperature,  and  density  in  addition  to  the  flow  velocit}-  over  the  channel 
surface  of  the  sample.  These  conditions  are  in  turn  used  in  an  equation 
that  e.xpresses  turbulent  heat  flux  to  a  flat  plate.  Ihc  heat  flux  is  com¬ 
puted  and  then  summed  to  yield  a  current  level  of  total  heat  input.  Ihc 
heat  flux  to  the  surface  is  also  applied  to  an  unstead>'  heat  conduction 
routine  in  which  the  surface  temperature  and  the  temperature  distribution 
in  the  test  sample  are  evaluated.  These  calculations  arc  all  performed 
within  the  same  time  step  of  the  overall  finite  difference  calculation. 
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The  technique  used  to  calculate  heat  flux  to  the  test  sample  surface 
requires  the  flow  velocity,  density,  and  viscosity.  The  density  and  visco¬ 
sity  are  evaluated  on  the  basis  of  a  reference  temperature  that  takes  into 
account  the  temperature  profile  resulting  from  the  boundary  layer  velocity 
distribution  and  the  sample  surface  temperature. 

Free  stream  velocity  is  defined  as 

U  =  [5  X  IO^’Yc  (T,  -  T  (A-18) 

where  T^,  the  free  stream  static  temperature  is 
Y-  1  2-1 

Tg  =  T^[l  +  C-^DNr] 


and  M,  the  Mach  number  of  flow  through  the  test  sample  is  given  by 


,Pl, 


y-  1 


M  =  Y 

Ps 


(A-19) 


where  equals  the  sonic  static  pressure  if  M  =  1,  or  the  downstream  pres¬ 
sure,  P2,  if  the  flow  is  subsonic. 


Density  and  viscosity  are,  respectively,  given  as 

-1 


“ref  '  Ps'SP,  • 


(A-20) 


P  =  ?  X  W-7  *  198)- 


1 


where  Tj.g£,  a  reference  temperature  is  defined  as 

-6  ^e" 

T  .  =  (T  +  T-..,„)/2  +  4.4  X  10  ° 
ref  e  SAMP  c^,Y 


(A-21) 


and  Tg^^ip  is  the  test  sample  surface  temperature. 

Turbulent  flat  plate  heat  flux  at  the  location  of  the  in-wall  thermo¬ 
couple  then  equals 


0.052(Pr^fUp«-%“-^  V  -  W 


(A-22) 


The  sample  surface  temperature,  TsAMp.  is  determined  from  the  one¬ 
dimensional  unsteady  state  heat  conduction  equation. 


3T  _  dh 
9t  "  “^3^2 


(A-23) 
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with  the  surface  boundary  condition. 


Q  -  .A-- 

for  X  =  0,  where  Q  is  the  heat  flux  to  the  sample  surface.  This  eciuation 
Joes  not  consider  the  effects  of  cylindrical  geometry. 

A  finite  difference  technique  using  a  geometrical  node  grid  spacing  was 
incorporated  into  the  Calspan  code  to  solve  the  unsteady  heat  conduction 
equation.  The  general  finite  difference  relationship  is  given  b>' 


(F  +  l)Ax-. 


(T. 


i-1 


—  T.  +  i.  T  ) 

F  i  F  i+r 


1-1 


(A- 25) 


where  o  is  the  thermal  diffusivity 
T  is  the  temperature  rise 
t  is  the  time  interval 
dxi  .  1  is  the  thickness  of  the  i-1  grid 

F  is  the  geometrical  multiplier  with  the  thickness  of  the  ith 
grid  being  F  times  that  of  the  i-1  grid. 


The  exposed  surface  boundary  condition  is  satisfied  b\’  first  determining  a 
fictitious  temperature  in  free  space. 


T  = 
0 


(A-26) 


This  temperature  is  then  used  to  establish  the  surface  temperature  rise  b\- 


2T^  +  TJ 


(A-2-) 


Barrel  Flow  and  Projectile  Motion 

The  flow  through  the  barrel  is  expressed  as  an  input  of  ma.ss  and  energy 
to  the  volume  between  the  test  sample  and  the  projectile.  As  mass  and 
energy  are  accumulated,  this  is  expressed  in  terms  of  pressure  and  tempera¬ 
ture,  which  in  turn  provides  the  accelerating  force  for  the  jirojcctile.  In 
this  calculation,  the  quasi-steady  assumption  of  the  previous  calculation 
is  relaxed  by  allowing  pressure  acting  on  the  base  of  the  projectile  to  be 
modified  according  to  the  Mach  number  of  flow  at  the  base  of  the  projectile. 
In  this  way,  the  unsteady  expansion  effects  of  flow  through  the  barrel  is 
taken  into  account.  The  equation  of  state  and  the  basic  energy  and  mass 
conservation  equations  are  the  same  as  those  for  the  test  gas  in  the  driven 
tube. 
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Projectile  motion  is  calculated  by  an  approximate  technique  that 
involves  determination  of  the  flow  Mach  number  at  the  projectile  base. 

M  =  Vj,[(Y  -  ljH2j'^/^ 


(A-28) 


where  Vp  is  the  projectile  velocity,  and  H2^,  the  static  enthalpy  at  the 
projectile  base,  is  defined  as 


Ht  =  2.5  X  10  H2 

^CO  ^ 


^  H2  - 


and  H2,  the  specific  enthalpy,  is  given  as 


where  Mg  is  the  mass  of  the  gas  contained  in  the  barrel. 

By  assuming  adiabatic  flow  and  ideal  gas  behavior,  the  isentropic 
relations  may  be  used  to  obtain  the  local  static  pressure  at  the  projectile 
base. 


Y  -  1  ,,?,-Y/Y  -  1 
^  =  P:„(l  +  — = —  M-) 


CA-29) 


where  P2  is  the  local  "reservoir"  pressure,  i.e.  the  static  pressure  of  the 
barrel  origin. 

Projectile  acceleration,  ap,  may  then  be  expressed  as 
a  =  52.17(p->  -  Pr)/i' 

where  Pj.  is  the  projectile's  local  resistance  to  motion,  and  Wp  is  the  pro¬ 
jectile  mass. 

Changes  in  projectile  velocity,  M'p,  and  displacement,  AXp,  can  be 
written  as 


AVp  =  apAt  and  AXp  =  VpAt  +  ap  — 


STG  Model  Printout 


(A-30) 


.A  sample  printout  of  the  STG  program  follows.  Definitions  of  "Initial 
Condition"  variables  are  given  in  the  program  printout.  However,  some 
explanation  of  the  "Tabulated  Output"  variables  is  required. 
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TIME 


PCH 


VOLl 


P1,T1,M1 

VOL2 


P2,T2,M2 

VELP,  XP 
VPROJ,  XPROJ 

WORK 


UE 

QFLUX 

QTOT 

TSAMP 


Time  after  piston  release  -  seconds 

Time  is  initially  zero  and  increases  to  a  maximum  value 
either  when  the  projectile  exits  the  barrel  or  when  it 
equals  a  limiting  value,  i.e.,  TFP, 

Driver  gas  pressure  -  MPa  (psi) 

PCH  is  a  prescribed  maximum  at  TIME  zero  and  decreases  as 
the  driver  gas  displaces  the  unlatched  piston. 

Test  gas  volume  included  from  piston  face  to  test  specimen 
inlet  -  m^(ft3) . 

VOLl  is  initially  the  entire  volume  of  the  driven  tube  but 
decreases  to  the  volume  of  the  test  gas  collection  chamber 
when  the  piston  has  been  fully  displaced. 

Test  gas  pressure,  temperature  and  mass  associated  with 
VOLl,  measured  in  MPa  (psi),  °K  (“R) ,  and  kg  (lb),  respec¬ 
tively  . 

Test  gas  volume  included  from  test  specimen  inlet  to 
projectile  base  -  m2(ft3). 

V0L2  has  a  minimum  value  of  the  test  specimen  bore  volume 
at  TIME  zero  and  increases  with  projectile  displacement 
to  include  the  entire  barrel  volume. 

Test  gas  pressure,  temperature  and  mass  associated  with 
V0L2,  measured  in  MPa  (psi),  ®K  (®R),  and  kg  (lb), 
respectively. 

Piston  velocity  and  displacement  -  m/sec  (ft/sec),  m  (ft). 

Projectile  velocity  and  displacement  -  m/sec  (ft/sec), 
m  (ft). 

Work  performed  by  driver  gas  on  the  driven  piston,  and  by 
test  gas  on  the  projectile.  During  piston  rebound,  the 
program  also  computed  negative  work  done  on  the  piston 
by  the  test  gas  -  J  (ft-lbf). 

Test  gas  free  stream  velocity  -  m/sec  (ft/sec). 

2  2 

Test  sample  surface  heat  flux  -  J/mm  -sec  (Btu/ft  -sec). 

Total  integrated  heat  input  to  test  sample  surface  -  J/mm- 
(Btu/ft 2) , 

Test  sample  surface  temperature  -  '’K(®R). 

TSAMP(l)  is  the  computed  temperature  on  the  test  sample  sur¬ 
face.  TSAMP(2)  through  TSAMP(9)  are  subsurface  temperatures 
computed  at  depths  printed  in  the  nonlinear  DELX  array,  which 
is  amended  when  the  time  increment  DELT  is  changed. 
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C 
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C  WPROJ  PROJECTILE  WEIGHT  -  LBM 
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C  TFO  LIMITING  TIME  -  SEC 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


74 


FORTRAN 

0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 

0011 

0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 

0021 

0022 

0023 


0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 

0055 


RELEASE  2,0  MAIN  DATE  »  81246  17/44/11 

READ  (5,101  )  NUMBER 
DO  99  NM=1, NUMBER 

DIMENSION  RG( 9)  ,CVG( 9  )  .CVTG( 9  )  .CVEXPGC  9  > .ALFAG( 9  > ,BETAG( 9  )  ,XG( 9  > 
DIMENSION  G1(9),G2<9),G3(9).G4( 9),G5( 9> 

DIMENSION  TNEUI 40  > ,TSAMP( 40 ) , DEL X<  40  ) , 0XNEW( 40  > . DTEMP <  40  ) 

REAL*4  M2,M20.M1 , KP  I 
00  =  1 
KK=1 
NN  =  0 

QFLUX=0.0 
UE-0.0 
R2=0.0 
CV2-0.0 
CVT2=0.0 
CVEXP  =  1 .0 
ALFA2=0.0 
8ETA2=0.0 
GAM2=1 .0 
P20=14.7 
T20=530.0 
XPRO0=0.0 
VELP-0.0 
VPR00=0.0 
C 

C  READ  IN  CONSTANTS  AND  VARIABLES. 

READ  (5,100)  PCH,VOLO,TO,CVOVR.RDVR 

READ  (5,100)  AT,KP1 .WSHELL.SHELRV.XPMAX 

READ  (5,100)  A2,POR12,C02,VOL1F,VOL20 

READ  (5,100)  BARL, BORED, PSTART.RESO.RESS.RESC.WPROJ 

READ  (5,100)  DELT.PRCI .TFO, FACTOR 

READ  (5,101)  NG 

WRITE  (6,103)  NM 

WRITE  (6,104) 

DO  3  ID=1,NG 

READ  (5,102)  XG( ID)  ,G1( I0),G2{ ID),G3( ID ) .(;4( ID ),GS( ID ) ,RG( ID ) ,CVG( 
♦I0),CVTG( IO),CVEXPG( ID>,ALFAG< ID>,B£TAG( ID) 

WRITE  (6,105)  XG(  I0),G1( ID),G2( 1D),G3( ID),64( ID),G5( ID) 

R2  =  R2  +  XG( ID  )*RG( ID) 

CV2  -  CV2  +  XG( IO)*CVG( ID) 

GAM2  =  GAM2  +  XG(  ID  )*( RG( ID  )/( 778 . ‘CVG! ID  )  ) ) 

ALFA2  -  ALFA2  +  XG( I D  )*ALFAG( 1 0  ) 

BETA2  =  BETA2  +  XG(  ID  )*BE TAG! ID  ) 

3  CONTINUE 

BORED  =  BORED/12. 

BORID  -  BOREO*12. 

ABORE  -  0. 785*BORED»*2 

CV=CVDVR 

DELTO-DELT 

DORI2  =  DOR12/12. 

DORID  -  DOR12»12. 

A12-0.78S4*DOR12»*2 
IPRC»IFIX(PRCI  ) 

IPRINT-0 

PCH=PCH*144. 

PCHP-PCH/144 . 

P20-P20*144 . 

PP-P20/144. 

PARTIM  •  0. 
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0056 

PSTART-PSTART*144. 

0057 

PSTARP-PSTART/144. 

0058 

QTOT-0.0 

0059 

RGAS-RDVR 

0060 

CP-RGAS/778.^CV 

0061 

6AM-CP/CV 

0062 

6M  -  PCH/<RGAS*TO)»VOLO 

0063 

E-GM*{CP-RGAS/778.  )*T0 

0064 

TIME  -  0. 

0065 

VMAX-XPMAX*AT 

0066 

XP-XPMAX 

c 

INITIAL  CONDITIONS  -  UPSTREAM  -  FROM  PISTON  FACE  TO  TEST  SAMPLE 

0067 

P1-P20 

0068 

V0L1-V0L1F+XP*A2 

0069 

T1-T20 

0070 

M1»P1*V0L1/R2/T1 

0071 

CVT2-0.0 

0072 

DO  1  J-l.NC 

0073 

CVT2-CVT2  +  );G<  J  )*CVTG<  J  )»(T1-460.0)»*CVEXP6(0  > 

0074 

1  CONTINUE 

0075 

E1-M1*<CV2'-CVT2)«T1 

0076 

RHOl-Ml/VOLl 

0077 

VVl-1 .O/RHOl 

c 

INITIAL  CONDITIONS  -  DOWNSTREAM  -  FROM  TEST  SAMPLE  TO  PROJECTILE  BASE 

0078 

P2-P20 

0079 

VOL2-VOL20+XPROJ*ABORE 

0080 

T2=T20 

0081 

M2-P2*V0L2/R2/T2 

0082 

CVT2-0.0 

0083 

DO  2  J-1  .NG 

0084 

CVT2-CVT2  +  XG(0  )*CVTG(0  )'»<T2-460.0)**CVEXPG<  0  ) 

0085 

2  CONTINUE 

0086 

E2"M2«(CV2*CVT2>*T2 

0087 

RH02-M2/V0L2 

0088 

VV2-1 ./RH02 

V. 

c 

PRINT  OUT  INITIAL  CONDITIONS. 

0089 

WRITE  (6,110) 

0090 

WRITE{6, 120)PP,PCHP,AT,T20,T0.KP1.R2,RDVR,WSHELL,CV2,CVDVR,SHELRV. 
♦CVT2.VOLO.XPMAX.CVEXP,PSTARP,ALFA2,BARL ,BETA2,A2,BORID,GAM2, 
♦OORIO,CD2,RESO,OELT.VOL1F.RESS,PRCI .VOL20,RESC,TFO,WPROJ 

c 

INITIALIZE  CONDITIONS  FOR  UNSTEADY  HEAT  CONDUCTION. 

009! 

DEPTH  -  0.5/12.0 

0092 

XK  -  19.3/3600. 

0393 

ALPHA  -  XK/57.6 

0094 

F-1 .3 

0095 

DELX0-SaRT(ALPHA»DELT/0.25) 

0096 

SUMX'DELXO 

0097 

0ELX( 1  )>OELXO 

0098 

DO  6  I*1 ,40 

0099 

TSAMP( I >-T0 

0100 

TNEW( I >-T0 

0101 

6  CONTINUE 

0102 

DO  7  1-1,39 

0103 

DELX(  I  +  l  )-DELX( I  )*F 
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0104 
0105 
0106 
0107 
0108 
0109 
0110 
01 1 1 
0112 
0113 
0114 
0115 
0116 
0117 


0118 

0119 

0120 

0121 

0122 

0123 

0124 

0125 

0126 

0127 

0128 

0129 

0130 

0131 

0132 

0133 


0134 

0135 

0136 

0137 

0138 

0139 

0140 

0141 

0142 

0143 

0144 

0145 

0146 

0147 

0148 

0149 

0150 

0151 

0152 

0153 

0154 

0155 

0156 

0157 


IFISUMX.GE. DEPTH)  GO  TO  7 
SUMX-SUMX  +  DELX(  I ) 

KTEMP  =  I 
7  CONTINUE 

I,ONSTl=ALPHA»2.0/(  1 .0  +  F  ) 

C0NST2"( 1 .0+F )/F 
C0NST3=1 .0/F 

TZIP-QFLUX“2.0*DELX< 1 )/XK+TSAMP( 2 ) 

10  CONTINUE 
T  =  T0 

RGAS-57 . 4+ . 333E-S»PCH 
CV=0. 175+. I83E-4*T 
CP=RGAS/77e.+CV 
GAM-CP/CV 
C 

C  COMPUTE  THE  DYNAMICS  AND  THERMODYNAMICS  OF  PISTON  MOTION. 
IF(VELP.LT.-0.1 )  WSKELL»SHELRV 
24  XO»XPMAX-XP 

DRAG=KP1*A^S< VELP ) 

FOVM-( ( PCH”AT-P 1»A2 )-DRAG)/WSHELL«32. 17 
26  DVELP  =  FOYM  •  OEUT 
28  VELP  =VELP  +DVEI.P 

IF IXP.LE. 0.00001 . AND. VELP. GT. 0.0)  V£LP=0.0 

DXP  =  0ELT*<  VELP  +  fOVM«0£LT/2.  1 

XP«XP-0XP 

IFIXP.LE.  .00001  )  GO  TO  99 
WORK  =  PCH  •  AT  •  DXP/778. 

E  «■  E  -  WORK 

30  VOL=VOLO+AT»XD 
RHOCH  =  GM/VOL 

T  =  E/{ (CP-RGAS/778.  )*6M) 

PCH  =  RHOCrl»RGAS*T 
C 

C  COMPUTE  TEST  GAS  PROPERTIES. 

CVTX=0.0 
CVTY-0.0 
DO  31  0  =  1, NG 

CVTX=CVTX  +  ( XG(0  I-CVTGIO  >*< (T2-460.  )»»CVEXPG(0  )  )  > 
CVTY=CVTY  +  ( XG(0  )*CVTG<0  >*<(Tl-460.  )*»CVEXPG(0  >  >  > 

31  CONTINUE 

CVY  =  ( (CV2  +  CVTV )»T1-ALFA2»RH01/778.0  >/Tl 
CVX»< ( CV2+:VTX  )»T2-ALFA2*RH02/778 . 0 )/T2 
RESIST  »  RESO  ♦  RESS  *  XPROO 

32  IF(RESIST.LT.RESC)  RESIST»RESC 
H2=E2*GAM2/M2 

H2INF»25000.*H2-VPROO**2/2. 

XM2INF»VPR0J/SQRT( (GAM2-1 .0)»H2INF  > 

FM-< 1 ,0  +  (GAM2-l .0 )/2.0*XM2INF**2  >••( -GAM2/( GAM2-1 .0  >  > 
DELP-P2»< 1 ,0-FM)/{ 1 .0+FM) 

33  P20-P2+0ELF 
P2INF=P2-DcLP 

PFORCE  -  A:,0RE  •  <  P2INF-RESIST) 

IFIPFORCE.LT.O.O)  PFORCE-0.0 

IF(P2.LT.PSTART)  PFORCE-0.0 

APROO  -  PFORCE/WPROJ»32. 17 

DXPROJ  ■  V7ROO  •  OELT  ♦  APROO  •  D£LT**2/2. 

35  XPROO  »  XPROO  ♦  OXPROO 

IFIXPROO.GE.BARL )  GO  TO  99 
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VPROO  -  VPROJ  ♦  APROO*  DELT 

V0L1-V0L1F+XP»A2 

VOL2-VOL20+XPROJ»ABORE 

36  PIAVE-Pl 
XMl-Ml 
XM2-M2 
IPASSaO 
P1SAV=P1 
P2SAV=P2 

37  CONTINUE 
IPASSaIPAS3+l 

PSTAR»P1*< 1 .0+(GAM2-l .0>/2.0>*»(-GAM2/<GAM2-l .01  > 

TTOT-Tl 

PS-PSTAR 

PTOT-Pl 

IF(P20.GT.P1 )  GO  TO  47 
IF(PSTAR.G;:.P20)  go  to  38 
PS-P20 

DELTFl1«8.02*A12*SQRT<  <P1-P20>*P1/R2/T1 >*D£LT 
GO  TO  49 

38  TSTAR^T1*2.0/<GAM2+I .0) 

0ELTM1=PSTAR'»A12*SQRT<GAM2*32.2/R2/TSTAR  >«DELT 
GO  TO  49 

47  CONTINUE 

PSTAR  =  PSTA  1*P20/'P1 
TTOT-T2 
PS=PSTAR 
PTOT-P20 

IF(PSTAR.GK.P1 I  GO  TO  48 
PS  =  P1 

0ELTMi=-8.  ./2*A12»SQRr<  <P20-P1  )*P20/R2/T2  )*DELT 
GO  TO  49 

48  TSTAR=T2»2 .0/(GAM2+l .0 ) 
0ELTMla-PSrAR*A!2»SQRT<GAM2*32.2/R2/TSTAR)*DELT 

49  CONTINUE 
XMlaMl-DEL VMl 
XM2»M2*DELTM1 
RH01=XM1/V0L1 
RH02=XM2/VaL2 
VV1=1 .O/RHOl 
VV2=1 . /RH02 

WORK--ABOR'  *OXPPvO«J*P2INF/778.0-P2/C33.*XPROO*3. 1416*BORED*DELT 
IF(VPROO,LV. 100. )  GO  TO  58 
C 

0ELH»CVY«DiLTMl'‘Tl»GAM2 

50  IF(0EI.TM1.LT.0.0)  OELH-  CVX*0ELTM1  »T2«GAM2 
EX»E2+WORKi-D£LH 

T2  =  EX/CVX/;!M2 

PX  =  (R2»T2  )/< VV2'SETA2)+ALFA2/VV2»*2 

P2=(PX+P2SAV)/2.0 

EY«El-DELri  +  PlAVi:/778.0*A2*DXP 

T1»EY/CVY/XM1 

PY-<R2*T1  ),'(  VVI-BETA2)+ALFA2/VV1**2 
Pl-(PY  +  P13..Y)/2.0 

54  DELPl-PY-P; 

PlAVE-Pl+OrLPl/2.0 
P20-P2  +  0EL.’ 

55  IF(  IPASS.lt. 3)  GO  TO  37 
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El-EV 

E2-EX 

P  1=PY 

P2-PX 

Ml-XMl 

M2=XM2 

XMOR-SQRT(  <  <  PTOT/PS  )•*<  <SAM2-1 . 0  >/6.^M2  >- 1 . 0  >  *2 . 0/ (  GAM2- 1  .  0  )  > 
TE=TTOT/< 1 .0  +  (GAM2-l . 0  ) /2 . 0«XMOR*«2 
UE  =  SORT(5.0t4*CYY*GAM2*(TTOT-T2  )  > 

TREF-<  TE+Ti;AMP<  1  >  )/2 . 0  +  4 . 4E-6»Uc  •*2.  CVY/r.AM2 
VIS=7.0E-7*TREF»*1 .5/<TR£F+198.0> 

RHOREF  =  PS/<  BETA>»PS  +  R2*TREF  ) 

QFLUX  =  0.05i»{RHOREF»UE )**0 . 8*VI S*»0 . 2»CVY»GAM2« { TTOT-TSAMP  (  1 ) ) 
IF(TSAMP{ 1 ) .GT. 3100.0 )  QFLUX»QFLUX»'  'TOT-SIOO. 0  ) / ( TTOT-TSAMP ( 1  )  ) 
QFLUX  =  QF'.UX'FACTOR 
T2IP  =  QFLUX-'2.0»DELX<  1  )/XK+TSAMP(  2  ) 

OTEMP< 1  )=A  .PHA/( DELX< 1 >**2 )•( T21 P -3 . 0«T3AMP ( 1 )+TSAMP(  2 ) > 

DO  56  K=Z,  HEMP 

OTEMPv  K  )=C  lilSTl  /  (  DELX(  K-1  )*»2  )*(  TSAM?{  K-1  ) -C0NST2*TSAMP  (  K  )  + 

•  CONST3*T3AMP(  ;c+l  )  ) 

56  CONTINUE 

DO  57  K=l.  TEMP 

TSAMP<  K  )  =  T3AMP( K )  +  0TcMP(  K )*DELT 

57  CONTINUE 

TSAMP(  KTEM  >  +  1 )  =  TSAMP < KTEMP- 1  > 

QTOT^QTOT+lr LUX-DEL T 
IF<  P20.GT.;M  )  U2  =  -UE 
GO  TO  60 

58  CONTINUE 
C 

V0I.  =  V0L1+V012 
XMT0T=M1+M'i 

CSUM=El+E2-Pl*A2*DXP/778.0  ♦WORK 
IU=XMTOT*v:jLl/VOL 
M2=XMT0T*V0L2/VJL 
E1=ESUM*V0T  i/VOi. 

E2  =  ESUM*V0.2/V0'. 

Tl  =  ESUM/CVi;/XMTOT 
T2  =  T1 

VVT=VOL/XMrOT 

?1>=<R2»T1  )/(  VVT-BETA2  >+ALFA2/VVT**2 
P2  =  P1 

RHOl-Ml/VO.'.  1 
RH02=M2/V0-2 
VV1=! .O/RHOl 
VV2=1 .0/RH02 
60  CONTINUE 

IF  <VOL1.GT.0.2^VMAX>  GO  TO  80 
IF( VOLl .LE .O.2*VMAX)OELT»O.l*DELT0 
IF(V0L1.LE.0.2*-/MAX.AND.KK.EQ.!)  GO  TO  62 
IF(  VOLl  .LE  .0.02  VMAX)DELT-=0.01»DELTO 
IF(  VOLl  .LE  0.02''VMAX  .ANO.KK.EQ.  2  )  GO  TO  62 
GO  TO  80 
C 

C  REVISE  GRID  S»ZE  DUE  TO  CHANGE  IN  TIME  STEP. 

62  0ELX0  =  SaRTvALr'HA>'DELT/0.25) 

SUMX-DELXO 
0XNEW<  1  >=d;;lxo 


FORTRAN  RELEASE  2.0 


MAIN 


DATE  -  81246 


17/44/U 


0269 

0270 

0271 

0272 

0273 

0274 

0275 

0276 

0277 

0278 

0279 

0280 

0231 

0282 

0283 

0284 

0285 

0286 

0237 

0288 

0289 

0290 

0291 

0292 

0293 

0294 

0295 

0296 

0297 

0298 

0299 

0300 

0301 

0302 

0303 


0304 

0305 

0306 

0307 

0308 

0309 

0310 

0311 

0312 

0313 

0314 

0315 

0316 

0317 

0318 

0319 

0320 

0321 

0322 


DO  64  1-2,40 

DXNEWl  I  )-DXNEW(  !-n*F 

IF( SUMX.GE .DEPTH )  GO  TO  64 

SUMX-SUMX+DXNEU( I ) 

KNMAX  -  I 
64  CONTINUE 

XKOLD-DELX< 1 ) 

XKNEW-0.0 

K»1 

KNEW-1 

A-(T2IP-TS.AMP(  2  )  )/<  2 . 0»0ELX(  lO  ) 

B-(TSAMP< 1  )-TSAMP<2)-A*DELX(IO)/DELX(X>*«2 
TNEW(  1  )-TS  ',MP(  1  ) 

66  KNEW-XNEW+1 
XKNEW-XKNEy+DXN£W( KNEW-1 ) 

TNEW<KNEW)=TSAMP( 1 ) -( A*XKNEW+B*XKNEV»*2 > 

IFIXKNEW.LT.XKO-DIGO  TO  66 

67  K-K+l 

XKOLO=XKOLO+DeLX(K) 

A  =  (TSAMP{ K-1  )-TSAMP( K+1  ) ) / ( ( F  + 1 . 0  )*DELX( K-1 >  > 

B  =  <TSAMP(K  )-TSAMP(K*l  )-A*DEUX{ K >  >/D£  .X( K  )**2 
63  KMEW=KNEW+. 

XKNEW=XKNEy  :-OXN:W<  KNEW-1  > 

XDEL-XKNEW'XKOLO+DELXIK) 

TNEW(KNEW)'  TSAM?(  K  )-( A*XDIL  +  B-XD£L*^-2> 

IF( XKNEW.L  .XKOlU .ANO.KNEW+1 .LE .KNMAa )  GO  TO  68 
J.F<  K+1  .  LE  .KTEMP  AND.  KNEW+1  .LE  .XNMAX  )  GO  TO  67 
DO  70  !=1,^NMAX 
DELX(  1  )-DX'.iEW<  1  > 

TSAMPC  I  )-T,(EW<  I  ) 

TNEW< I )=T0 
70  CONTINUE 
KTEMP-KNMA/ 

KK  <=  (CK+1 
80  CONTINUE 
C 

C  PRINT  OUT  COM'>UTED  RESULTS. 

IF  (00. EO. •  )  GO  TO  88 

87  NN=1 

PARTIM=PARriM+DELT 
IF( TIME .GT.TFO)  GO  TO  99 
IF(  IPRINT.LT.  IPilC  )  GO  TO  90 
GO  TO  89 
C 

88  WRITE  (6, IB7) 

WRITE  (6,103) 

00=00+1 

89  PCHP-PCH/1..4. 

P1P=P1/144.0 
P2P-P2/ 144  . 

TIME-TIME  +  i’ARTIM 
PARTIM  =  0. 

WRITE! 6, 190  >  TIME , P IP .T1 .Ml , VOL  1 .VELP , VPROO , WORK , QFLUX , PCHP . P2P , T2 
♦,M2,VOL2,XF ,XPROJ,UE,QTOT.(TSAMP( I >.I»1 ,9) 

1F( NN.EQ.O  )  GO  TO  87 
IPRINT  «  0 

90  IPRINT-lPR.NT+1 
GO  TO  10 
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0329 

0330 

0331 


0332 

0333 


0334 


99  CONTINUE 

100  FORMAT  (8F10.5) 

101  FORMAT  < 1 1 > 

102  FORMAT  ( F5 . 2 , 5A3 , 6E 1 0 . 3  ) 

103  FORMAT  (12H1  INPUT  DATA . T48 TEST  NO.',TS9.in 

104  FORMAT  { / / //T6 , *  TEST  GAS  MIXTURE  BV  MOLE  FRACTION'/) 

105  FORMAT  ( /T6 , F6 . 3 . T1 7 , SA3 ) 

110  FORMAT  < ///T6, ' INITIAL  CONDITIONS  FOR  TEST') 

120  FORMAT  <//T6,'P20  =>  •  ,T1 7  .  F  1 3 . 6 ,  T3  i  .  '  PS  I  •  .T48.'PCH  =  '  ,  T59 .  F  1 3 . 6 . 
+T73,'PSI'  ,T90,‘AT  =  '  ,  TlOl  .  F  1 3 . 6 .  T1 1 5  ,  ‘  F  T>**2  '  / /T6  .  '  T20  ='.T17. 
♦F13.6.T31 . 'DEGR* .T48. 'TO  •’ ,T59 , F 1 3 . 6 , T73 . 'DEGR' .T90. 'XPI  » ' ,T101 , 
♦F13.6.T115, 'LBF/FT'//T6. ' R2  » ‘ . T1 7 . F 1 3 . 6 . T3 1 . ' FT-L EF /L BM-DEGR ' . 
♦T48, 'RDVR  = ‘ ,T59.F13. 6,T73, ‘ FT-L BF /L 3M- DEGR ' .T90.'WSHELL  ='.T101, 
♦F13.6.T115. ■LBMV/Te. •CV2  »*.T17.F13.6.T31, ‘ BTU/LBM-DEGR ' ,T43. 

+  'CV0VR  «• .T59.F13.6.T73. ‘ BTU/L BM-DEGR ' .T90, 'SHELRV  « '  , T  1 0 1 . F 1 3 . 6  . 

♦  T115, 'LBM'//T6, 'CVT2  = ' , T1 7, F 1 3 . 6 , T48 . ' VOL 0  =  '  . T59 . F 1 3 . 6 . T73 . 
♦'FT**3' ,T90, 'XPMAX  = ' .TIO 1 , F 13 . 6 . T 1 1 5 . ' F T ' //T6 , ' CVE XP  =' ,T17. 
♦F13.6.T48. ‘PSTAR7  = ' .T59 . F 13 . 6 . T73 . ' P S I '  //T6.'ALFA2  = '  . T 1 7 . F I  3 . 6 . 
+T90,'BARL  =■ ,T10! ,F13.6.T115. 'FT' //T6. 'BETA2  = ' . T 1 7 . F 1 3 . 6 . T48 . 

♦ ‘A2  = ' .TS9.F13.6.T73. •FT**2' .T90, ‘BORED  ='.T101.F13.6.TI15.'IN’// 

+  T6,‘GAM2  =■ ,T17,F13.6.T48. 'DORIZ  » ' . T59 . F 1 3 . 6 . T73  .  '  1 N  '  / 

♦  /T48.  'C02  =•  ,T59,F13.6.T90,  'RESO  =  '  .  T 1 0 1  ,  F  1 3 . 6  ,  T 1  1 5  .  '  L  B-l '  / /T6  . 

+  'DELT  =■ .T17,F13.6.T31 . 'SEC* .T48. 'VOL  IF  = ' . T59 . F 1 3 . 6 , T7 3 . ' FT**3 '  . 

♦  T90. 'RESS  ='.T101.F13.6.T115, ' LBM/FT ' //T6 . 'PRC  I  = ' , T 1 7 . F 1 3 . 6 , T4  8 . 
♦' VOL20  =' .T59,F 13.6.T73, 'FT**3' ,T90. 'RESC  = ' . T 1 0 1 , F 1 3 . 6 . T 1 1 5 . ' L BM ' 

♦  //T6, 'TFO  =' ,  ri7,F13.6.T31 , 'SEC' .T90. 'WPROO  =  '  . T 1 0 1 , F 1 3 . 6 . T 11 5 . 
♦'LBM'/////) 

187  FORMAT  ( 18H1  TABULATED  OUTPUT) 

188  FORMAT! ///T6. 'TIME ' ,T20, 'PI ' .T34. ' T 1 ' . T48 . ' M 1 ' . T62 . ' VOL  1 ' , T76 . ' VE L 
+  P ' ,T90, 'VPROO ' ,T104, 'WORK' . T1 1 8 , ' QFl UX ' / 16 . ' PCH ' . T20 . ' P 2 '  . T34  .  ' T2  ' 
♦,T48, 'M2' . T62, ' V0L2' .T76. 'XP '.790,  ’ XPROJ ' . T1 04 ,  '  UE  '  .  T 1 1 3 .  'QTOT'/ 

♦  T4. 'TSAMP< 1 > ‘ .T18. ‘TSAMP<2) ' ,T32, 'TSAMPl 3  > ' .746,  'TSAMP( 4 ) ' , 

♦  T60. 'TSAMP( 5  ) ' .T74 . ' TSAMP( 6  > ’ .T88. 'TSAMPf  7 ) ' .T102.  'TSAMP( 8  > ' . 
♦T116, 'TSAMP(9)‘// ) 

190  FORMAT! 9E 1 4 . 6/9E 1 4 . 6 /9E 1 4 . 6/ ) 


0335 

0336 


C 


STOP 

END 


TEST  GAS  MIXTURE  BY  MOLE  FRACTION 


TABULATED  OUTPUT 


^  o  o. 
u.  wx 
ao< 


oom  OOU1  ooin 


ooot  oou)  oou)  oom 


OOO  900  ooo  ooo 


oom 

ooo 


OuiX 

3  3  < 


O  O 
O  O  -- 
ei  oc  a. 
a.  a.  Z 

>  X  < 


(A 


uj  o.  Z 
>  X  < 


oom  oom  oom  oom 


>  o 

o 

ooo 

ooo 

o 

o  o 

o  o 

o 

OOO 

o  o 

o 

ooo 

ooo 

OOO 

o 

o 

o 

O 

o 

© 

o 

o 

c 

o 

o 

o 

PD 

PD 

(D 

t*l 

rs 

rs 

PD 

CD 

PD 

rs 

CS 

PS 

rs 

O 

o 

o 

o 

o 

o 

o 

o 

o 

O 

— 

o 

-- 

♦ 

♦ 

♦ 

♦ 

♦ 

♦ 

* 

♦ 

4- 

4 

4 

4 

Ul 

Ul 

Ul 

Ul 

Ui 

UJ 

Ul 

Ui 

Ui 

Ui 

1 1 1 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

m 

m 

m 

m 

m 

m 

m 

m 

m 

m 

m 

PD 

CD 

rs 

PD 

(n 

PD 

PD 

PD 

PS 

PD 

rs 

PS 

rs 

i  O 

m 

oom 

oom 

o 

o  m 

o  o 

m 

oom 

o  o 

m 

©  o  m 

0 

0 

m 

oom 

oom 

o 

o 

LO 

O 

o 

m 

o 

o 

m 

>  O 

o 

ooo 

ooo 

o 

0 

0 

o  o 

o 

ooo 

o 

o 

o 

ooo 

o 

O 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

•  N 

PD 

N  (N  (D 

pD  <N  rs 

PS 

(N  <n 

rs  (N 

PD 

PD  rg  PD 

(D 

(N 

PD 

PD  rg  rs 

PD 

(N 

rs 

cs 

(VJ 

PS 

rs 

PN 

rs 

n 

(N 

rs 

(S 

rs 

rs 

)  O 

o 

OOO 

ooo 

o 

o  o 

o  o 

O 

ooo 

o 

o 

o 

ooo 

O 

O 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

•  ♦ 

♦ 

*  *  * 

♦  ♦  ♦ 

♦ 

*  * 

♦  *■ 

* 

♦  ♦  ♦ 

♦ 

♦ 

♦ 

♦  ♦ 

4- 

4 

4 

♦ 

♦ 

♦ 

4 

4 

♦ 

4 

4 

4 

- 

cu  w 
0»  ^•k  o 
m  o 
moo 
^  o  m 
m  o  («) 
mom 


^  ^  o 
O  fv  o 
moo 
^4  m  m 
o  o  <«i 
morn 


o  M  o 
O  (SI  o 
rs  QD  o 
(n  so  m 
-« <n  <*» 

>*  rs.  m 


o  o 
moo 
o 

m  (sj  m 
m  o  (n 
—  fs.  m 


00  rs  o 
m  (V(  m 
m  (n 
(SI  m 


^  o 

(Si  (n  o 

(>i  (SI  o 
(•>  CD  m 
m  (^  (n 
(SI  m 


«Si  (»>  o 
fs.  o 
—  —  o 
os  o  m 
o  O  (*) 
CO  f»»  m 


m(M  o 
(SI  o  o 
o  m  o 
(SI  m  m 
m  fs*  CO 
<*)  m  m 


m  CD  o 
sD  m  o 
os  o 

(SI  N.  m 
os  c^  n 
(n  m  m 


Os  sO  LfS 
N  0\  (»1 
^  m  m 


rs  o>  o 
m  rs  o 
0^0 


m  m 
o>  o  n 
»  m  m 


m  r*.  o 
O'  r»  o 
m  c*.  o 


css  m  o 
^  rs  o 
m  (>.  o 
o  r*.  js 
m  O'  rs 
m  M  m 


ooo  ooo 


ooo  ooo 


ooo 


ooo  ooo 


ooo 


ooo 


o  o  z 
>  >< 
m 


fSI  o  o 
ooo 
c-s  O  O 
oom 

^  N  rs 
(M  CD  m 


(SI  o  o 

(S|  O  O 

«.  o  o 
(B  o  m 

^  csi  rs 

(N  CD  m 


♦  I  * 

UJ  Ul  Ul 

moo 
ooo 
(SI  o  o 
mom 
^  (SI  rs 
(SI  00  m 


(SI  o  o 
1^0  0 
(BOO 

oom 
^  (N  cn 
(M  CO  m 


ooo 
•"  o  o 

CO  o  o 
•»  o  m 

rs  (SI  rs 

(SI  00  m 


ooo 
o  o 
rs  o  o 
o  l(S 
ca  (SI  rs 
(SI  00  m 


moo 
r-  o  o 
(^  o  o 
CD  o  m 

—  (SI  CO 

(SI  CO  m 


moo 
0S  o  o 
ooo 
CD  o  m 
o  (SI  rs 
(Vi  00  m 


ooo  ooo 


U(  Ui  IW 
CD  O  O 
O  O 

moo 
mom 
O'  (SI  rs 
-•  tio  m 


(Si  o  o 
N  o  o 
CO  o  o 
mom 
CO  r;  r» 

-•  CD  m 


—  o  © 
ooo 
ooo 
O  O  iis 
rs  rj  r> 

—  CO  m 


rs  o  o 
r*  o  © 

—  o  o 
mom 
m  rj  rs 

—  CO  m 


moo 
(SI  o  o 

-9  0  0 
^  ©  m 
rs  (s,  n 
-•  CO  m 


moo 
^  o  o 
OD  O  O 
N  o  m 
rj  ri  r* 
—  CO  US 


ooo  ooo 


ooo  ooo 


ooo 


ooo  ooo  ooo  ooo  ooo  ooo  ooo 


ooo  ooo 


(SI  Z 

z  z  < 


♦  I  ♦ 

W  Ul  Ul 

m  (SI  o 
(SI  (Vi  o 
o  «•  o 
CB  rs  m 

—  (»S 

IN  fs  m 


♦  I  ♦ 
Ul  Ul  Ul 

—  o  o 
(N  (N  o 
o  m  o 
CO  o  m 

—  <N  (^ 
(SI  m 


♦  1  ♦ 

Ui  Ul  Ul 

(n  «  o 

—  (NO 

ooo 

oom 

—  (N  « 

(N  rs  m 


o  rs  o 
o  m  o 

O  (N  O 

PS  (N  m 

—  rs 

CN  rs  m 


Ui  Ul  Ul 

m  m  o 
m  m  o 
os  M  o 
rs  m  m 
—  o  rs 
(N  rs  m 


rs  (N  o 
<N  o  o 
os  rs  o 
rs  o  m 


♦  I  ♦ 

W  Ul  w 

rs  rs  o 
CO  rs  o 
(o  m  o 
rs  o  m 
—  m  cn 

(N  CD  m 


^  (N  O 

rs  o  o 
CO  o  o 
rs  o  m 
—  o  (^ 
Psi  os  m 


->•  rs  o 
rs  V  o 
rs  ^  o 
rs  —  m 
—  rs  rs 
Pj  os  us 


^  o 
ooo 
moo 
fs  m  m 
—  o  rs 


N  • 


rs  m  o 
os  CO  o 
moo 
rs  ^  m 


'  m 


UJ  Ul  Uj 

rj  m  O 

moo 

CO  O 
ps  ts  m 
—  (SI  ('S 
N  —  US 


OOO  ooo  ooo  ooo 


CS  (N  C 
CD  ^  O 

PJ  o  o 
US  m 


ooo 


ooo  ooo 


ooo  ooo  ooo  ooo 


•  (N  Z 

•  h-  < 


^■9  0 

m  m  o 
VO  m  o 
oom 

(*>  CD  CO 

m  m  m 


♦  ♦  ♦ 
UJ  Ul  UJ 

^  o 

rs  (s  o 
PS  PS  o 

—  —  m 

rs  CO  cs 

m  m  m 


m  m  o 
«  ^  o 
m  m  o 


<N  (N  o 
09  CO  o 

os  os  o 
oo  CD  m 
rs  po  rs 
m  m  m 


CO  CO  o 
m  m  o 

--NO 

m  m  m 
^  ^  (*s 
m  m  m 

ooo 


(n  (n  (*> 
ooo 
♦  ♦  ♦ 
Ui  Ul  Ui 

CO  ro  o 
O'  c^  o 
M  o 
<n  cn  m 
m  m  cn 
m  m  m 


Ul  Ul  Ul 

rs  fs  o 
m  m  o 
(N  (N  o 
<n  c*s  m 

VO  so  (*S 

m  m  m 


♦  ♦  ♦ 

Ul  UJ  Ul 

m  'O  o 
oo  CO  o 
m  m  o 
m  m  m 
PS  rs  rs 
m  m  m 


Ul  UJ  Ul 

-s  —  o 

ooo 
m  m  o 
oom 
os  (7>  r> 
m  m  m 


oo  CD  m 
o  o  rs 
m  m  m 


CO  CO  o 
(N  (N  ■T' 

ooo 
oom 
rs  P'S  CO 
m  m  m 


m  m  o 
o  O  'I 
—  o 
mom 
m  m  rs 
m  m  m 


O'  CTS  o 
rs.  ps  m 
50  CD  rs 
m  m  m 


cr 

rs  rs  o 
rs  fs  vO 
(N  CN  rs 
rs  rs  in 


ooo 


ooo 


ooo  ooo 


ooo  ooo 


ooo 


ooo  ooo 


ooo  ooo 


ooo 


.  a.  < 

m 

V- 


—  o 
OS  os  o 
rs  rs  o 
ps  rs  m 


^  ^  o 

ooo 

CO  CD  O 

(0  CO  m 

^  ^  rs 

—  —  m 

ooo 


ooo 
«« o 
rs  ro  o 

—  —  m 
m  m  (n 

—  —  m 


m  m  o 
CO  00  o 
rs  n  o 
m  m  m 
m  m  cs 
—  —  m 


CO  o  o 
CO  CO  o 
o 

—  —  m 

VO  m  n 

—  —  m 


m  m  o 
os  o>  o 
CO  CO  m 
m  m  rs 
—  —  m 


(N  (N  (n 

ooo 
♦  ♦  ♦• 
Ul  UJ  Ul 

n  rs  o 
ps  Ps  o 
ooo 
OS  OS  m 
rs  rs  n 
—  —  m 


UJ  UJ  UJ 

CO  <B  O 
eg  (N  o 
ooo 
(N  (N  m 
os  os  po 
—  —  m 


ooo 

K  ps  O 

m  m  o 

CD  CO  m 
O  O 
(N  (N  m 


♦  4-  ♦ 
Ul  Ul  Ui 

m  m  o 
rs  rs  o 
p^  rs  o 
ov  Os  m 
eg  M  rs 
N  (N  m 


M  M  o 

m  m  o 

(N  (N  O 
rs  ps  m 

m  m  rs 
eg  (N  m 


(N  (V  o 
rs  rs  o 
^  ^  o 

(^  (*>  m 
os  O'  cn 
pg  (N  m 


^  —  o 

—  —  o 
^  ^  o 
eg  eg  m 
^  m 

rs  PO  m 

ooo 


ooo 
«  —  o 
CO  »  o 
oom 


ooo 


ooo  ooo  ooo  ooo  ooo 


ooo 


ooo 


ooo  ooo 


Z  Z  a. 

—  u  z 

a.  < 


('sw  — enrs  — rsrs 

oo  ooo  ooo 

♦  ♦  I  ♦  ♦  1  ♦  ♦ 


—  (v»  rs 

ooo 


—  rs  rs  -•  rs  m 
ooo  ooo 


UJ  UJ  Ul 

(NO  o 
o  o  o 

—  o  o 

Ps  m  o 

CB  rs  — 

o  m  m  — 


UJ  UJ  Ui  Ul 

o  o  o  rs 

mo  o  IN 
rs  o  o  (N 

^  m  o  rs 

os  PS  CD 

mm  (N  m 


Ul  Ul  Ul  UJ 

O  O  (B  o 

o  0^0 

O  0  —  0 

m  o  m  m 

PS  -•  ps  (n 

m  rs  m  m 


W  Ul  Ul  UJ  Ui  u 

CJS  rs  O  O  ps  O 
O'  PS  O  O  O'  o 
OS  VO  O  O'  (N  O 

O'  —  m  os  ^  m 
o  m  (n  o  ^  PO 
^  m  m  m  m  m 


Ui  Ui  Ui 

Os  ^  o 
o*  ts  o 
os  fs  o 
Os  ^  m 

O  (N  (^ 

m  m  m 


—  PS  PS 

ooo 

I  *  * 
Ul  Ul  UJ 

O'  ps  o 
ooo 
OS  fs  O 
o*  (s  us 
O  O  P'S 

ps  m  m 


—  PS  rs 
ooo 
1  ♦  ♦ 
Ui  Ul  UJ 
CD  IS  O 
OS  CO  o 
OS  m  o 
O'  —  m 
o  CO  rs 
CO  ^  m 


—  PS  P'S 

ooo 

I  *  * 

Ul  Ui  UJ 

rs  m  o 
os  P'S  o 
c^  pg  o 
O'  O'  m 
o  m  n 

O'  ^  vT) 


o  m  PS 

ooo 


O  (N  o 
o  os  o 
o  rs  o 

—  mm 
o  P'S  rs 

—  ^  m 


o  rs  PS  o 
ooo  o 


O  rs 
O  O 


Ui  ul  Ul  Ul  VU  Ul 

ooo  0—0  om 

0  —  0  omo  ossf 

orso  o(NO  osm 

—  ^m  — r'm  opj 

—  —  rs  tN<3'PS  rsfs 

-.^us  — rsm  — 


Ooo  ooo 


ooo  ooo 


ooo  ooo 


8.^ 


.J 


3000 


140999E*00  0.510935E*02  0.777261E*03  0.2170S2E*01  0 . 1  .•S524E « 02  0.573502E»03 


r 


84 


'*S9t*00  0.420254E^03  0.l449y5E*04  0 . 2  1  1  4  I OE -^0  I  0.2:/094E*01  0.614660E*03 


^  ^  C^f 

ooo  ooo  ooo 


\0^^O  CDO'u'J 

(sj^o  ^4^i)o  n>«o  o-ro  ^ 

omo  lOCDo  omo  ni~>o 

incBiA  MiAtfl  (Njtsij) 

^  (>i  cD«<*i  flocoin 

—  tn  nr*.i/)  ^criji 


^  r\j  m  w  ^  n  yr  f*> 

COC  ooo  OOi^  ooo 

♦  ♦♦♦  ♦♦♦ 
UJ  ^  4^  4*J  4^  4^  ^ 

COOO  ooo  yf  CC 

fMJio  ®ors.  fsjof**- 

ocsjo  «co«  ooco 

ijinm  ocorn  i/inn  c*io\n 

ooin  <voi/>  j?  —  tO  ®  — 1/> 


ooo  ooo 


ooo  ooo  ooo  ooo 


fo  —  r)  r>  —  M  m— icn  ^.-fn 

ooo  ooo  ooo  ooo 

+  *•  ♦♦♦  ♦♦♦  *  ^  * 


<*)ino  0^0  ONO 

Ovom  r>»oo  or>>OD  co«<si 

^400  ^ooo 

ro*«u)  «co« 

4O«0C*>  i/>  — 1*»  niflr) 

lOtsiiA  lAcsjin  i/)^tn 


(MOO  ^  fv 

»  —  irt  ®  n  — 

rs.  r*.  (s.  0  05^ 

0(^0  O  O  10 

—  (MV  0—10 

0  —  0  0-0 


ooo  ooo  ooo  ooo 


rvon  — <M0  foof^  (MOO 

o  o  o  o  p*»  r>  o  ®  ®  CP  o 

—  0(n  (MO0  —  —  O  (M  —  o 

(V,  M  T  «■  -.  T)  1^—0  •«»  O 

^  <•  •  •  ••  •  —  - • 

OO0  f>»OD0  0O0  VO^> 


o  V  —  ®  n  o 

0  SO  0  O  0  0 

0  V  0  O  ®  (M 

-  T>  0 


ooo  ooo  ooo 


OOO  ooo 


I  b84E 


94959t»00  0.5ZObo2E«04  0  ..'90844E +04  0.20139JE+01  0.J.4935E+00  0.493821E+03  0.7508b9E‘03  -  0  .  1  6  064  4E  -  0  1  0.612647E+04 


*  ♦ 

M  UJ 
O' 
7'  xO 
Si  r> 
Si  ^ 
N 

-I  .0 


3  O 


•*4  LU 
O  CO 
D  Ji 
NJ  <Sl 


S»  (SJ 
SJ  ^ 
a  CO 


S'  o 

\i  iT) 
S)  kSi 
■n  ^ 
3  r*- 

- 


s  o 


X 

T*  O' 

a  N 


J'  N. 

o  o 


*S. 
30  rs 
—  vO 
r>  fn 


O'  m 
^  "V* 
cn  j'. 


^  ^4  M 

o  o  o 
♦  ♦  ♦ 


^  Ui  Ui 
CO  fO  fVi 
xO  O'  xo 
X/)  CD  O' 


^  xO  lil 
—  rg  W 
o  xn  xo 
O'  03  m 
rv  m  in 
xO  — 


O  O'  CD 
« cn  in 
in  00  (SJ 
m  (SJ  K 
in  xn  in 
rs.  «  xn 


i  O  CD 
>  xn  m 
)  —  n 


>  CO  Oi 
r  xo  in 
5  m 


«•  xO  rs 

in  in 

—  (•»  PS» 
O'  rxx  o 

O'  K  xO 
00  —  xn 


in  -* 
in  CD  xD 
O'  —  xn 


o  o  o 


o  o  o 


o  o  o  o  o  o 


s0  «  O 
(SI  O'  C'' 
xO  CD  Ifi 
xO  (SJ 

o  O'  n 

—  —  in 


o  o  o 


GO  (*) 

o  o  xn 

(SJ  m 


fn  CO  in 
M  m  o 
in  in  (*) 
in  CD  in 

M  M  in 

—  (M  in 


ro  o  xn 
(SJ  n  n 

—  (SJ  ^ 


O'  —  e 
^  rn  O' 
O  «  (SJ 
(Si  (*)  CO 


N  M  — 

—  ^  JSJ 

w  cn 

—  rv  O' 

«  xn  xo 

—  (SJ  n 


o  o  o 


o  o  o 


o  o  o  o  o  o 


o  o  o  o  o  o 


o  u*) 

(SJ  O' 
(Si  rv 
xn  m 
xo  <n  rv 
xsj  in 


Ul  UJ  ^  w  i 


n  • 


O' 


00  C+J  fv 
O'  CO  PI 
O'  xn  CO 
xn  xn  rs 
—  (SI  n 


xn  (SJ  <■» 
in  p-  (SJ 

O'  (SI  o 

O*  CO 

fs.  xn  so 

—  (SI  in 


CO  (SI  ^ 
xn  psi 
M  CD  p*> 

o  o  n 

CO  CO 
M  (M  in 


-•  O'  ^ 

xn  w 
xn  « xn 
xn  m  xn 
CO  rs>  CO 
-•  (SJ  xn 


—  xn  (Tx 
O'  p'.  or 

—  (M  in 


r*.  O'  (sj 
O'  K  O' 
—  (SI  m 


o*  —  — 

(SJ  xD 
M  in 
^  (Si  in 
O  CO  O' 
(SI  (SI  in 


—  pn  CD 
U)  so 
M  CD  ox 
(SI  pg  in 


(S'  ^  (Sj 

■n  —  — 
xn  n  (SI 
CO  CD  Psi 
—  CD  O 
(SI  (SJ  in 


O'  O  M 
CO  xn  in 
rs.  o 
xO  pg  in 
(Si  O'  o 
psj  (sg  xn 


—  n  (». 

M  in  in 
xO  xO  « 

xn  in  O' 
P^  O'  o 
(Si  (SI  vn 


xn  O'  po 
^  O'  — 

(SJ  rsi  xD 


o  o  o 


n  (SJ 

xn  O'  o> 

P>.  O'  — 

M  If 

3  0  0 


O  O  O 


o  r*»  pn 
O'  o  Pg 
•n  rg  in 


o  o  o 


rv  ^  ro. 
xn  in  CD 
O  —  (SI 
CO  (SI  xO 


o  o  o 


cn  (SJ  « 

M  m  PI 

CO  (SJ  fv 

o  (*>  (n 

(SI  (SJ  fO 

00  (SJ  gs 


o  o  o  o  o  o 


ooo  ooo  ooo  ooo 


o  o  o  coo 


Ox  m  ps 
CO  (*3  in 

iD  00 

n  cn  cn 

CD  pj  xn 


O  —  O' 
{'I  fs.  o 

(SI  O' 


«  O'  xn 
CO  m  o 
in  in  — 
<n  00  o 


(n  (SI  crx 
xn  o  <n 
o  m  o 


CD  o  n 

O'  o  — 
xn  PJ  (^ 
N.  in  (s. 
O'  xo  xn 
CD  PJ  xn 


in  in  ^ 
—  «n  p*+ 

xn  PJ  CO 

(n  m  00 

fv  lO 
O'  PJ  xO 


ooo 


O'  ^  in 
fsj  CO  r** 
O'  N  xn 


ooo 


PJ  lO  ^ 
•  (SJ 
CO  o  o 
xn  PI  m 
^  O'  CO 

os  PJ  xn 


ooo  ooo 


in  xn  (SJ 
xn  CO  p- 
PI  (n  -n 
PJ  m  o 
xn  o  O' 
O'  po  xn 


UJ  X«J  UJ 
^  PI  pn 
V  O'  in 
pg  o  rv 

P^  —  O' 

O'  (*)  xn 


CO  O'  o 
^  O'  r*. 
ox  in 

OO  Pv  — 

»  CO  CO 
^  (Tx  -n 


P«.  —  <SJ 
fs.  — .  (Sj 
fs.  in  O' 
cn  CO  CO 
CO  M  x'i 

^  O'  n 


rv.  —  n 
ps  n  pg 
(SJ  xn  sO 
CO  o  n 
(s.  O'  o» 
^  CO  n 


ooo 


(SJ  M  ^ 

p>.  o 

^  CO  p'. 


ooo  ooo 


xn  —  — 
xO  —  O' 
—  (SI  — 
xo  xn  PI 
xn  O'  — 


xn  M  o 

pH  CP  — 
^  O'  ^ 
O'  O'  PJ 
in  <M 


!S»  9*  — 

xn  o  cn 

^  pH  (>. 


vO  PJ 
O  O'  O' 

If  If  CO 
^  O'  — 

^  m  ^ 

^  xD  pH 


ooo 


ro  o  po 
ooo 
♦  ♦  ♦ 
Ui  UJ  1>J 
O'  pH 
O  (SI  PO 
(O  pH  CO 
If  m  ^ 
n  —  lO 
^  xn  pH 


ooo 


fO  O  PO 

ooo 
♦  ♦  ♦ 
Ui  Ui  Lkj 
pH  O  O' 

»  (O  pH 

po  in  CO 
ph  (SJ  (SI 
PJ  pH  in 
^  in  pH 


ooo  ooo 


—  O  (SI 
in  pH.  fs. 
xn  «  in 

pH  O 


(SJ  (St  pH 

—  PO  O' 
pH  oi  in 
O  CO  CD 


ooo 


—  «  in 
CD  O' 
in  (SI  •" 
in  O'  o 
O'  ^  o 

P3  *»  CO 


ooo 


n  o  ^ 

CO  —  — 
«  ^  CO 


ooo 


m  Ox  o 

'*H  M  — 

»  m  m 
•»  —  ■(« 


ooo 


x^  Id  UI 

xn  ^  r» 
lO  O' 


PI  O'  pH 
OOO 


m  O'  •h. 

ooo 


ooo 


^  PI  O' 
PJ  xo  n 
(SJ  xn  xn 
■n  ps  X 
n  3  -H 
PI  O'  pH 


•h.  m  rs. 

o  xn 
(SI  rn 

^  pH 

O  Ph  <»5 
PJ  —  CO 


—  m  O' 

—  m  — 
o  xn 

PJ  -*  03 


OOO 


Ph  OC 
(M  PJ  in 

o  CO  PI 

—  xf  pH 

O  Ph  sn 
PJ  —  CO 


Ph  in  in 

in  ^  ^ 

j3  Ph  ^ 
O  XT'  C 

PI  O'  ^ 

OOO 


OOO  ooo 


Ph  in  — 

PH  ^  — 

MOO 


ooo  OOO 


Pi;  —  -  —  .Tx  • 


-« O'  m 
P3  O'  in 

—  CO  P-» 

•X  O 


in  ^  xn 
O'  xn  xo 
P3  m  (SJ 


UI  UJ  Ui 

in  o  CD 

O  ^  O' 
O'  PI  o 


ooo 


P3  O'  O* 
CO  (SI  m 
If  O'  O' 

Ph  (SJ  p. 


ooo  ooo 


CO  O'  O' 
^  (SI  Ph 
Ph  if  CO 
«  pp 


«  CO  CD 

*•  so  pH 

CO  o  ^ 
O'  J“  -T 


M  O'  » 

ooo 


PI  —  » 

ooo 


M  —  O 

ooo 


(M  —  CD 

ooo 


(SJ  —  O' 

ooo 


—  J3  O 
P3  O  Ph 

o*  ®  PS 

hh  'O  Ph 
PI  —  O 
OOO 


o  « (n 

O  pH  ffl 
M  —  « 


O'  PJ  O' 
^  O'  in 
Ph  xn  ® 
o  os  o 
o  pH  o 

M  —  O' 


CO  ®  XJ 
O'  O'  s 
xn  o  - 
o  —  c 
o  ®  • 

(SJ  —  c 


J  m  ps 
»  O'  PS 
P  xo  \D 
S  (SI  CO 
5  o  cn 
J  —  O' 


®  ' 


r  PJ 


xf  O'  pH 
PJ  PS  PS 
O  ^  O' 

o  ®  in 

PJ  —  O' 


UJ  Ui  Ui 

xO  •»  © 
OOO 

o  (SJ  in 
o  ®  — 
o  ®  ® 

PJ  —  O' 


}  o  ^ 
>  PJ  o 
'  o  (n 
'  O'  o 


>  PS  O' 

)  PJ  us 

t  ^  PH. 
'  M  US 
'  O'  O 


^  O'  — • 
N  O  pH 
(SI  o  ® 


f  —  — 

xn  o  o 

CO  pH  — 

0'  O'  — 


ooo  ooo 


ooo  ooo 


O'  US 

—  ® 


Ui  Ui  Ui 

SH  C^  -t 

O'  O  •■•J 

-H  —  n 

us  m  a 
O'  o  — 

Pi  M  O' 


OOO 
♦  ♦  ♦ 
UJ  Ui  Ui 
us  O  — 


O'  ®  xf 

a  o  PS 
M  PJ  a 


O  Ph  a 

o  'js  a 
us  PS  ^ 
^  Ph  xf 

o  o  xn 
PI  PJ  a 


us  a  xo 
rn  us  ^ 
^  ®  PS 

a  a  pH 

O  O  Ph 

PS  PJ  a 


ooo 


»  a  PJ 
®  o  xn 
^  o  a 
—  —  a 
PS  (SJ  a 


ooo 


ooo  ooo 


OOO  ooo 


a  O  Ph 

■w  a  — 

—  in  PS 


uJ  Id  UJ 
—  (SO 
a  cn 
a  ^  a 

us  PS  ^ 
PJ  —  o 
PS  PJ  — 


ooo 


Id  Ui  Id 

am® 
n  a  — 

o  —  ^ 

m  PJ  a 


•H.  rs  — 

a  PJ  o 

M  rn  rn. 


UI  ■ 


r  a 


OOO 


—  in  a 

o  in  a 
®  us  o 
o  a  ^ 
PJ  PJ  o 


ooo 


in  • 


ooo 


pH  »  a 
O  ^  PJ 
®  m  a 
^  us  a 
o  (SI  o 


ooo  ooo 


a  PJ  -• 
ooo 


♦  ♦  ♦ 
Id  Id  Ui 

—  pH  ® 

PJ  —  PJ 
PJ  m  ^ 

o  pH  o 


a  PJ  — 
ooo 


®  (SI  ^ 

If  a  M 
o  a  us 
^  a  a 


a  (SJ  — 
ooo 


♦  ♦  + 
Id  Id  UJ 

a  m  — 
.n  M  JS 
a  xn  a 

—  PJ  ® 

in  PJ  — 
a  M  — 


^  ^ 
ooo 
♦  ♦  ♦ 


a  xn  o 
xn  (SJ  pg 
a  pg  — 


a  m  a 
—  a  r 
a  ®  pg 

f  00 

®  PJ  pg 

a  (SJ  — 


pg  ' 


ooo  ooo 


ooo  ooo 

♦  ♦ 


PJ  pH  pg 

a  ®  (SI 

m  PJ  — 


^  PJ  pH 

a  — 
o  o  a 
Ph  a  us 
O  PJ  — 
®  —  — 


♦  ♦ 

I  Ui  UJ 

a  pH 
PJ  a  a 
Ph  us  ® 
®  a  a 

®  PJ 


a  • 


OOO 
♦  ♦  + 
xd  Ui  UJ 

^  pH  — 

a  a  o 

—  (SJ  PJ 

a  —  (SJ 
a  a  PJ 
a  —  — 


a  —  f 
pH  PJ  c 
o  a  xj 
oaf 


®  xn  PJ 

PJ  Ph  ® 

o  a  — 
a  xn  o 
o  a  a 


«  o 

(SI  —  — 

®  pH  V 

—  a  a 


ooo  ooo 


ooo  ooo  ooo  ooo 


xn  a  a 
o  a  a 
JS  xn  n 
a  ®  o 


'  PJ  PJ 

1  a  a 
■  M  a 
s  m  a 
'  ®  o 


•  (SJ  - 


♦  ♦  ♦ 

Xd  Ui  Ui 

®  m  US 
xn  a  xn 
a  o  ® 
us  m  ^ 
a  ®  — 


Ui  Xd  Ui 

®  a  o 
in  a  a 
^  a  ® 


.  (\j . 


o  a  ^ 
ooo 
♦  ♦  ♦ 
Id  UJ  Xd 

o  us  o 
xn  o  PJ 
us  a  M 
xn  ^  — 
a  xo  PJ 


•  PJ  - 


o  a  ^ 
ooo 
♦  ♦  ♦ 
Ui  Ui  Ui 

®  a  PJ 
xn  xn 

Ph  Ph  P^ 
US  Sf  ® 

a  ®  Pg 


.  pj  . 


®  o  xn 
m  a  xo 
CO  m  a 

us  T  PJ 

a  ®  a 


CO  —  a 
in  ^  a 
a  us  ^ 


Xd  Xd  Xd 

®  xn  ® 
us  m  pH 
o  ^  a 


Xd  id  Xd 
®  « 
in  pH  US 
—  a  ® 


®  a  xn 
xn  ®  xn 
(SJ  (SJ  a 


-  PJ . 


ooo  ooo 


ooo  ooo 


iS^S 


M  in  — 
®  a  a 
xn  ®  ® 


ooo  ooo 


—  (SI  pH 

pH  us  o 
pH  a  »» 
in  a  pH 
u^  o  — 
PJ  a  ® 


pH  pH  — 

®  o  pg 
(Si  a  ® 


o  —  a 
o  Ph  a 
us  a  (SI 

®  (SJ  Ph 

a  PJ  o 
a  a  pH 


ooo 


a  pH  (SI 

PJ  m  a 
mao 
a  PJ  a 

®  pH  (SI 

a  a  o 


ooo 


a  Ph  o 
(SJ  m  a 
^  m  ® 
•*  0  o 


—  —  a 

ooo 


CO  —  o 
®  ®  ® 
in  —  — 

®  o  m 
a  o  — 

—  N  — 


® 


M  X  r 

O  PJ  ' 

us  PJ  ® 

Ph  m  pg 

a  PJ  — 


ooo 


m  PJ  a 
®  a  a 

.a  PJ  — 


ooo 


ooo 
*  *  * 
xd  Xd  Xd 
O  Ph  a 

®  m 

pH  m  ^ 

®  o  a 

PJ  ^  a  9*  ♦  ^ 


a  pH  ® 
a  m  a 
®  us 
o  a  ’ 


>  m  a 
»  pH  o 


>  a  ® 
s  o  — 
s  iM  xn 


98386£«01  0.ie5358E*00  0.354547E»03  0.)014)9E*04  -0 . 280203E -0 1  0.162517E»05 


oo  ooo  ooo  ooo  ooo  ooo  ooo  ooo  ooo  OOO  OOO  OOO  OOO  OOO  ooo  ooo 


Ui  Ui 
N  O 
a*  CD 
so  — 

M  \0 

—  <VJ 

r> 


^  o  at 
o 

at  —  M 
fsi  n  \o 


(n  ^ 

(M  ^  ^ 
(*)  to  at 
OOO) 
-•  <sj  r» 
M  n  «o 


CD  c*t  to 
to  M  O' 
O  CM  — 
<D  lA 
CM  CM  ^ 
m  n  to 


1  *  * 
Whl  kiJ  UJ 

CM  rv  at 
^  at 
at  00 
ho  ot  to 

••  CM  ^ 
fO  C*»  to 


—  ^  rt 

ooo 

I  ♦  ♦ 
Ui  Ui  IM 

to 

—  cn  ot 
CM  to  so 
O  V  CM 
r**  n  m 
cn  c*i  to 


lO  (M  to 
c*)  to  to 
o  at 
tO  ot  CD 
at  m 
cn  cn  to 


o  cn  w 

CM  (M 
(Tt  cn  «o 
Pa  cn  to 
CM  ^  to 
^  cn  lo 


C*>  f*.  CM 
tO  ^ 

^  cn  to 
ooo 


«  OD  at 
«  o  — 

tO  tO  CM 

P^  o  o 

00  to  CD 
«  C^  to 


to  to  (*> 

CO  P^  CM 
^  O  CM 

«  c»t  a> 

CM  lO  CD 
lO  ro  lO 


to  CD 
lO  ^  cn 
CM  CM  p.. 

n  ^  to 
to  to  at 

lO  P^  tO 


o 

^  ot  a» 

ro  pv 
lO  tO 
O  to  o 
sO  ro  p** 


CM  CO  CM 
lO  —  Ni 
M  to  tD 
®  fn  tO 
®  to  (Si 
O  M  P^ 


O'  — 

pg  01  ( 
p.  n  f- 


O  O 


OOO 


OOO  ooo  ooo  ooo 


ooo 


OOO 


OOO 


OOO 


ooo 


ooo 


pg  tn 
O  CO 
CO  to 
CO  — 
CO  p^i 

o  o 


*•  *  * 
LU  UJ  LU 
P>.  »0  CM 
CO  vO  to 
CM  pg  at 

CO  CO  i/J 
O  ^  CM 
«*  CO  p». 


pg  CM  CM 

«  Ps  P*. 

—  to  ^ 
in  CO  to 
o  tn  CO 

—  CO  P^ 


CM  Ot  to 
m  CO 
O  CM  O 


^  —  fO 
OOO 
♦  +  ♦ 
ui  Ui  ua 
sO  ^  CO 
tO  CO  ^ 

O  O  to 
at  tn  CD 
o  p<»  in 

—  CO  p». 


^  ^  o 

to  at  to 
^  o  at 
^  to  o 

««  00  CM 
—  CO  P^ 


tn  at  CM 
in  -sp 

PI  PI  — 
PI  p>.  ^ 

^  at  CD 
—  CO  p«. 


CM  CO 

tn  to  CM 

«a  p>‘  CM 
tn  CD  CD 
—  O  <T» 

^  P«» 


in  PI  CD 

to  in  tn 
o  ^  PJ 

CD  O  CO 
—  CM  ^ 
M  «  CD 


« tn  CD 
O  CD  O 
tO  CO  po 

o  CM  at 
CM  CO  pg 

M  «  CD 


tO  P»  P* 
Pw  p^  ^ 

pg  in  9 

P)  «  tO 

PJ  ♦  ♦ 

-a  «  O 

ooo 


CO  ^  — 

O  o  to 
to  P*  ^ 
N  in  to 

«  O 


pg  Ot  to  CO  to  at 


‘  at  o 

O  P.  O' 
CM  »0  CO 


CO  tn  to 

—  CD  pg 

pg  pj  « 
CO  00  o 
^  ^  ot 


in  at  — 
tn  CM  ^ 
CO  CM  tO 
tn  o  j) 
CO  O'  pg 
—  tr  O' 


tc  o  a 
CO  O'  P» 
CO  O  M 

—  in  at 


ooo 


ooo  ooo 


ooo 


ooo 


ooo 


ooo  ooo 


ooo 


ooo  ooo 


68B42E 


127l50Et03  U.MZI40E<04  -0 . 769636E -0 1  0.Z04173E«0S 


I 

I 


*  * 
i4i  UJ 
00  -• 
o  o 

lO  (Tk 

o  a 
VO  rg 

>0  vO 


U)  «  <o 
*<V  ot  vO 
^  09  m 
kO  <l>  iO 
CD  r»  (M 
VO  so 


vO  n 
O  O  O 
*  ♦  ♦ 


r*.  ^l*  o 

vO  ov  cn 

vO  VO 


u»  M  m 
e  o  o 
♦  ♦  ♦ 
UJ  Ui  M 
U)  4}  rN 

vn  K  o 

CO  -M  rv 

o»  m  ui 


*  *  * 
Ui  UJ 
VO  U) 

M  r*  o^ 
me-f" 
<n 

M  (VJ  w 
«-  r*»  vo 


lO  N  n 
o  o  o 


eO  M  fs. 

m  vo  u) 

lO  CO  o> 

<rt  so 

—  M  ^ 

M  hs  so 


^  sn  vn 
K  m  os 
VO  O  OS 
so  (vg 
O  U)  lO 
r>.  vo 


••  N  fn 
o  o  o 
♦  >  ♦ 
W  Ui  M 
^  ^  ^ 
VA  (Si  O 
09  -N  (SI 
(Si  •-  o* 
vA  sO  VA 
c^  r>.  sO 


^  O  ^ 

V  en  <o 
vA  o  ee 

(Si  O  tA 
lA  rs»  so 
a  (s..  >0 


«  eo  (Si 

sO 

(So  rs  tO 


Ui  ui  ui 

(Ss  n  o 
—  (*»  fs.^ 
(S.  iA  OS 
GO  oo  rv 
so  r»  sO 


«  (Si  m 

e  e  o 


CD  (SI  fl» 
00  N  (*) 
(SI  00  OS 


(*1  ^g  O 
K  GO  vO 
GO  vO  (*) 
OS  (<>>  « 
(A  OS  OS 
j)  rs  vo 


rs.  SO  isi 

M  OO  « 
VA  OS  OS 
rv  (SI  — 

o  e  o 

lA  GO 


—  M  n 
os  O  so 

(*9  GO  VO 
M  N.  Os 
SO  O  O 
«  (B  rv 


so  (S  « 
(SI  OO  OS 
vO  ^ 
O  (Si 
(M  — 
«  00 


o  o 


ooo  ooo  ooo  ooo  ooo  ooo 


o  o  o 


ooo 


ooo  ooo  ooo  ooo  ooo 


ooo  ooo 


so  (Si 

(SI  -* 
OS  CO 
N.  rs. 

*« 

o  o 


^  M 

ooo 

I  ♦  ♦ 
Ui  Ui  Ui 

o  e  GO 

(A  VA  OS 
vA  00  (SJ 
(^  OS  (SI 
^  (A  rs. 
CO  (*)  K 


^ 

(SI  o  (SI 
OS  (*9  (n 
VO  SA  lA 

r*.  (n  ® 
o  (n 


—  ^  m 
ooo 

I  ♦  ♦ 


rs.  OS  o 
VO  vA  o 
rs  o  ® 
o  (n  OS 

OS  m  rs 


I  ♦  ♦ 
Ui  W  Ui 

o  ®  ^ 

rs  cn  VO 

®  ®  VO 

sA  m  (Si 
n  (Si  — 

os  (n  ® 


—  os  rs 

00  o  n 
vA  ® 
M  eO 
vO  (Si  (SI 
OS  Cl  o 


I  ♦  > 
Ui  Ui  Ui 
Cvi  eA  ® 
•M  «o  (n 
o  «  eg 
^  vO 
GO  ^ 
©  09  ® 


(A  OS  rs 
O  M  lA 
O  iA 

(n  00 


lA  VO 
OS  vO  OS 
^  lA  (^ 
(SI  IS  O 
O  O  IS 

—  m  « 


(*)  ds  00 
Os  w  IS 
-»  C9  OS 


«  OS  rs 
O  OS  «o 
K  M  ^ 
tA  OS  OS 
O  OS  OS 
M  (SI  GO 


Ui  Ui  tM  Ui  k 


o  ^  rs 
rs  ea  p9 
e  os  M 
—  (Sj  os 


r  rs  —  4 


OS  ®  -• 
O  ® 

-•  (Si  Ov 


M  rs  «■ 
O  ^  vA 

—  «  lA 

—  N  O' 


OOO 


ooo 


ooo 


ooo 


ooo  ooo  ooo 


ooo  ooo 


ooo 


ooo  oeo  ooo 


o  o 

N.  rs 
OS  OS 

Oi  o 

(SI  rs 
VA  OS 


ooo  ooo  ooo  ooo  ooo  ooo 

*  *  *  *  *  *  ^  ^ 

Ui  Ui  Ui 
rs  VA  vfl 
m  VA  o 
o  rs  M 
^  VA  ^ 


N  O  vfi 
N.  eO  O 
VO  C9  vA 
VA  ®  ^ 

cn  ov 

VA  ov 


O'  (^  •S 
^  VA  O 
•H  VA  -• 


eO  VA  OS 
(M  -s  fs 
OS  (SI  (M 
(M  O  ^ 
vA  VO  O 
VA  — 


VA  (SI  OS 
•«  C^  ® 
vO  sO  vO 
VO  (*9  vO 
tA  ®  O 
M  VA 


®  « 
N  (*1  VO 
m  vA  o 
o  ov  c^ 

vO  OS  O 
—  vA  — 


^  VO 


OOO 
♦  ♦  ♦ 
Ui  Ui  Ui 

VA  ^  vA 
^  ^  vA 

rs  tn  vo 
rs  (SI  <n 

VO  (*)  M 
M  so  ^ 


vo  n  Os 
(A  O  (Si 
^  (n  CO 
-«  os  VA 


ooo  oeo  ooo  ooo 


<•«  o  vo 
O  (*)  o 

M  EO  C^ 
VA  lO  rs 
rs  VO 

•m  ^ 


(SI  M  ® 
(SI  rs 

rs  (*9  CO 
O  «  os 
rs  a  — 
mm  iO  ^ 


M  ^  rs 

(SJ  IS  (»9 

(n  n  rs 
(SI  (SI  •« 

a  o  (Si 

M  rs  «« 


lA  (^  O 

a  a  a 
a  IS  « 
A  O  (*9 

a  (Si  (Si 

•M  rs  M 


ooo 
*  *  * 
Ui  VM  UI 
rs  « rs 

A  «  o 

(^  A  — 
a  a  A 
a  C9  (SJ 

-« rs  M 


ooo 
*  *  * 
VU  Ui  UI 
C9  —  rs 
a  A  n 
rs  eO  A 
(Si  a  A 
a  A  (Si 

M  rs 


A  a  » 
—  o  o 
A  a  a 

a  rs  |Vg 

M  rs  — 


o  o 


ooo 


ooo  ooo 


ooo 


ooo 


ooo  ooo  ooo 


ooo 


ooo 


ooo  ooo 


ooo  ooo 


A  A 
A  o 
a  — 

rs  Os 
O  (SI 


o  rs  ov 
rv  (Si  A 
«  A  A 
a  (SI  (SI 

A  (SI  (^ 


vu  Ui  Vii 
rs  OS 
rs  «■ 
(SJ  n  — 
a  o  A 


a  < 


’  C9 


A  o  a 

—  A  <*1 
N  — 


*  *  • 
Ui  VU  Ui 
^  (SJ  A 

«  o  a 

rs  A  (*5 

(9  a  A 

A  O  V 

N  (Si  — 


A  C9  <9 

a  ••  (SI 
O  (SI  o 
^  fs  a 
a  (*)  ^ 
(SJ  eg  — 


A  rs  ^ 
(Si  A  O 
O  A  A 
(9  (SI  « 


<9  0^ 
OOO 
♦  ♦  ♦ 
Ui  Ui  Ui 

A  a  o 

O  O  A 

eg  o  A 
a  a  cg 
M  a  A 

(9  (Si  *• 


rs  o  eg 

4V  A  A 

(9  a  « 

^  o  ^ 

(9  (9  A 
(9  (9  — 


a  A  A 
sv  A  A 
(9  (9  — 


A  A  O 
A  a  GO 
AAA 
O  O  fs 
A  O  A 
(9  ^ 


(Si  a  (9 
(9  A  eg 

—  «r  a 

(SI  (S  CD 
IS  pi^  A 
(9  ^ 


o  a  A 

^  ©  rs 
A  (9  a 
N  A  a 
a  rs  A 

a  ♦ 


a  rs  A 

a  eg  rs 

—  a  ^ 

—  (9  — 
O  A  A 

«  A 


OOO  OOO  OOO  ooo  OOO  OOO 


ooo  OOO 


ooo 


VU  VU  Ui 

—  a  eg 
eg  (9  ^ 

-■  a  (SJ 
m  —  s, 

—  f 


(SI  <9  a 

A  a  <9 
a  (SJ  — 


lul  UI  V«i 

a  A  ^ 

A  A 
O  O  — 


—  <9  A 

a  o  o 
(9  eg  ^ 

—  Ip  ® 

n  f.i  s. 


^  O  eg 
a  a  « 

^  <9  (SI 

a  >0  — 


A  A  eg 
©  a  <9 
(9  A  A 

eO  rs  f-« 


o  o  ^ 
ooo 


rs  —  a 
rs  a  A 
^  (9  A 


rs  a  es 

n  rs  A 
a  ©  — 

(-1  O  es 
A  (^  » 


O  O 

ooo 


a  a  a 

(Si  A  A 


A  A  (SI 
O  a  (9 
A  A  eg 
(9  (SI  a 
®  (9  X 


0  0  9 

ooo 
♦  ♦  ♦ 

VU  Ui  ViJ 

a  •-  o 

—  (9  eg 
o  o  a 


o  o  « 
ooo 


O  O  «V 

ooo 


(SJ  ^  (9 
*-  O  rs 

eg  a  N 

A  VO  O 


o  o  « 

OOO 


A  O  ^ 
eg  (9  © 
ANN 
rs  a  o 
f.  -  a 


^  (9  A 
^  a  ts 
rs  A  a 
a  a  a 


M  ■ 


A  (Si  (9 

W  A  — 

a  —  A 
—  —  a 


ooo  ooo  ooo  ooo  ooo  ooo  ooo  ooo  ooo  ooo  ooo  ooo  ooo 


—  ^  A 
a  A  -v 
A 

a  —  A 
a  a  o 

—  eg  (SI 


rs  A  o 
a  a  fs 
a  a  a 
a  A  a 
a  a  o 
—  eg  (SJ 


a  (SJ  (SI 
a  o  -4 
-•  (9  eg 


a  a  9* 
«  A  A 

a  ^  — 
rs  rs  A 
a  o  -4 

-4  (9  (SI 


—  O  4» 

ooo 
♦  ♦  * 
ViJ  UI  ViJ 
A  —  A 
a  o  a 

4»  —  rs 
rs  (SI  A 
a  —  — 

—  (9  eg 


♦  ♦  ♦ 
Ui  VU  iU 

a  A  fs 

A  A  (SJ 
o  (9  a 
rs  A  fs 

a  —  — 

—  r9  (SJ 


a  A  (SI 

AAA 
A  (SJ  <9 

A  o  a 
a  eg  -4 

—  <9  (SI 


O  (9  A 

—  (9  <9 
<9  a  sr 

A  (9  a 
a  (SI  — 

—  <9  (SI 


— •  A  a 
a  »•  4« 
a  -4 

A  rs  ® 

a  eg 
—  (9  eg 


a  A  A 
rs  (9  A 
A  »  A 
A  O  es 

a  (9  -• 

—  (9  (M 


-4  C  -V 

OOO 

♦  ♦  ♦ 


^  a  rs 
A  eg  A 
a  (9  — 
—  (9  eg 


A  a  a 
o  A  a 
a  a  ^ 

^  rs  » 

a  (9  — 

—  <9  (Si 


A  —  A 
a  A  IS 

A  O  >• 
4V  O  (9 

a  ^  — 

—  (9  (SI 
OOO 


—  o  'e 

ooo 
♦  ♦  ♦ 
lii  VAi  A 

a  —  V 

fs  M  rs 
«  —  rs 
^  eg  — 
0  ♦  — 

—  (9  (Si 

ooo 


Vii  A  Ui 

a  eg  A 
a  N  o 
eg  o  (9 

»  O 

a  ^  — 

-4  (9  eg 


ooo 


ooo  ooo 


ooo 


ooo 


ooo  ooo 


ooo 


ooo  ooo 


ooo  ooo 


ooo 


ooo  ooo 


O  A  A 
—  in  a 
a  (9  (9 
a  eg  — 

A  A  ^ 
^  (9  eg 


—  a 
-4  rs  (9 
a  A  ^ 

4»  (9  (SJ 


-4  a  (9 
(Si  o  A 
A  fs  ^ 
^  (9  (Si 


V  mt  ^ 

ooo 
*  *  * 
UJ  lU  Vii 

o  a  (SI 

o  (9  ^ 
(M  (9  A 


♦  ♦  ♦ 
A  UI  Vli 

^  (Si  a 
A  a  eg 

>«  —  a 

(SJ  A  A 
<9  fs 
(9  (M 


+  ♦  ♦ 
A  A  VU 

pg  —  a 
a  —  o 
a  a  -4 
—  rs  a 

pg  rs  ^ 
^  a  eg 


ooo  ooo 
♦  ♦  ♦  ♦  ♦  ♦ 


A  —  ^ 

.N  a  o 
a  a  a 
—  a  a 


a  a  — 
o  —  eg 
—  A  a 
(SI  a  eg 
O  rs  ^ 

^  a  eg 


AO©  ooo 
♦  ♦  ♦  ♦  ♦  ♦ 

U  W  U  A  U  U 
ej  rs  fp 

—  A  A 
A  es  ^ 

eg  a  o 

a  es  ^ 
a  a  eg 


a  o  A 

a  o  rs 
^  rs  a 

a  rs  a 
a  a  eg 


♦  ♦  ♦ 
tii  VAi  Ui 

a  a  eg 
«  A  A 
(s  a  — 
^  a  A 
ps  fs  a 
a  a  eg 


^  a  a 
A  a  a 
A  rs  a 
a  a  (SJ 


ooo  ooo 


o  rs  a 

A  rs  (SI 
a  a  eg 


rs  —  vo 

a  A  A 

^  IS  (Vg 

a  a  eg 


A  Vii  Vii 

A  (SI  a 
rs  vT 
a  o  (SI 
A  A  « 
a  ps  (SI 
a  a  rs* 


ooo  ooo 


ooo  ooo 


ooo 


ooo 


ooo 


ooo 


ooo  ooo 


^  a  o 

A  a  rs 

a  (Vi  pg 


♦  ♦  * 
VU  tu  Ui 
»  pg  -A 
a  A  A 
rs  ei  ^ 
a  A  4^ 
(SI  a  ps 
a  eg  rsi 


AAA 

a  a  A 
a  *4  o 


a  "s 

A  A  a 
rs  o  A 
pg  a  isg 


^  ^  A 

rs  o  o 
(M  A  a 

o  o  rs 
A  *4  A 

pg  a  eg 


♦  ♦  ♦ 
vu  A  vu 
a  -4 
A  A  >4 
a  (SI  a 
A  a  A 
eg  —  A 

eg  a  pg 


♦  ♦  * 
Ui  Vii  vu 

A  pg  a 
o  (SI  ^ 

A  O 

A  A  eg 
o  >4  a 
pg  a  eg 


A  eg  (vj 
A  a  a 
a  a  a 
A  A  a 
a  -4  A 
-4  a  (SI 


♦  ♦  ♦ 
VU  Ui  vu 

a  pg  o 
a  Ps  a 
o  (Si  o 
a  a  a 
a  4*  a 

—  a  (SI 


♦  ♦  ♦ 

Ui  lu  VAi 

^  eg  a 
A  A  o 
ps  a  es 
a  A  eg 
A  «4  A 
—  a  CM 


ooo  ooo 


a  44  eg 
-4  «  a 
o  A  a 

^  -4  ^ 

—  a  (M 
ooo 


o  —  a 
.i4  a  a 
>4  a  a 
a  a  A 

(M  -4  ^ 

-4  a  psi 


Ps  pg  pg 

•M  M  ^ 

-4  a  (M 

ooo 


rs  o  a 
o  -4  a 
44  a  pg 


^  a  a 
^  a  (SI 
A  a  pg 


a  rs  a 

—  a  -9 

—  a  — 
a  A  a 
44  o  a 
a  a  pg 


ooo  ooo  ooo  ooo 


ooo  ooo 


ooo 


ooo 


ooo 


ooo 


ooo 


ooo  ooo 


ooo 


rs  a  A 
A  rs,  (M 
-4  A  a 


A  O  44 

Pi  a  rs 
a  a  eg 


A  4*  a 

a  a  ^ 
a  a  o 
a  A  a 

-4  (M  (M 


♦  ♦  ♦ 
vu  VAi  vii 
rs  a  ^ 
A  a  a 
4»  A  a 
a  a  fs 

a  A  a 

44  eg  eg 


IS  a  o 
A  (Sg  A 
A  PS  a 
a  a  a 
a  a  a 

44  pg  (M 


fs  a  4e 
A  rs  o 
a  rs  a 
a  a  a 

a  a  rs 
—  eg  (SJ 


+  4  ♦ 

VaI  Ui  VaJ 
fs  a  4f 
A  a  eg 
rs  a  a 
a  a  A 

a  a  rs 
-•  eg  pg 


4  4  4 

A  VAi  VaI 

rs  eg  a 
A  a  a 
a  a  a 
a  a  o 
a  A  rs 
—  (M  eg 


4  4  4 
AAA 
rs  ^  ^ 

A  A  a 
a  a  A 
a  a  a 
a  a  a 

—  <Sl  (M 


4  4  4 
AAA 

Ps  a  A 

A  44  Ps 

0  0  4* 

a  ^  (M 
a  A  A 
—  eg  (M 


AAA 

Ps  a  rs 
A  a  o 
44  o  a 


ps  a  a 

AAA 

(M  44  rs 


rs  <M  44 
A  a  • 
a  (M  rs 

a  ^  a 
a  a  ^ 
—  eg  eg 


rs  a  eg 

A  O  Ps 

^  a  a 


rs  a  a 
A  a  A 
A  a  — 

a  ^  (M 

a  A  4* 

—  (SI  pg 


A  a  a 

A  «V  A 

a  ^  a 

—  pg  (SI 


ooo  ooo  ooo  ooo 


ooo 


ooo  ooo 


ooo  ooo 


ooo  ooo 


ooo  ooo 


ooo  ooo 


S8 


i 


•JU 


09E^0]  0.2*:4406£^00  -0.409404E 


72E 


APPEN'DIX  B 


ISENTROPIC  EQUILIBRIUM  COMBUSTION  CODE 


91 


FHSCBOUn  BUIK-JW®  llli® 


APPENDIX  B 

ISENTROPIC  EQUILIBRIUM  COMBUSTION  CODE 


Adiabatic  Compression  Program 

An  existing  equilibrium  combustion  program  has  been  modified  to  calculate 
the  properties  of  an  STG  test  gas  mixture  at  various  points  along  an  isen- 
trope.  The  computer  code  uses  an  internal  listing  of  constituents  and  their 
thermal  properties  to  dictate  the  instantaneous  composition  and  thermal 
properties  of  the  STG  test  gas  at  any  specified  pressure. 

Input 

Input  consists  of  an  initial  temperature  and  pressure,  along  uith  25 
other  pressure  values.  The  pressure  schedule  ranges  from  atmospheric  to  500 
MPa  in  20  MPa  steps  but  the  input  pressures  must  be  specified  in  atmosjiheres . 
Also,  for  every  constituent  in  the  STG  test  mixture,  a  "reactant"  card  must 
be  submitted  that  gives  the  program  its  initial  composition  and  mole  frac¬ 
tion.  "Omit"  cards  may  also  be  submitted  which  exclude  the  indicated  species 
from  consideration  in  any  reaction  processes.  Solid  carbon,  or  C(s)  is 
routinely  omitted  from  STG  runs  because  of  the  high  C(s)  concentration  alread) 
present  in  the  specimen  surface  which  would  inhibit  any  further  formation. 

Operation 

The  initial  mixture  and  mixture  ratios  are  analysed  to  determine  which 
chemical  elements  are  present.  The  number  of  gram  atoms  of  each  initial 
constituent  is  calculated  along  with  the  total  molecular  weight  and  cnth:ilp>’ 
of  the  mixture. 

The  taped  thermal  data  is  searched  and  the  names  and  thermal  properties 
of  possible  compounds  which  could  be  formed  from  the  available  atoms  are 
extracted.  This  list  is  then  compared  to  the  "omit"  cards  and  the  net  list 
of  species  to  be  considered  is  printed. 

The  initial  composition  is  then  varied  in  an  attempt  to  minimize  the 
Gibbs  free  energy  of  the  mixture.  Species  with  mole  fractions  of  less  than 
10"7  at  any  iteration,  are  dropped  from  consideration.  The  initial  numiier 
of  atoms  of  each  element  present  must,  of  course,  remain  the  same. 

The  program  calculates  the  entropy  of  the  initial  composition,  winch 
is  then  held  constant.  For  the  next  value  in  the  pressure  sclicdulc,  a 
corresponding  temperature  is  estimated.  Based  on  this  estimated  temperature 
a  new  equilibrium  composition  and  mixture  entropy  arc  calculated.  If  the 
new  enthropy  and  the  fixed  enthropy  agree,  the  program  advances  to  the  next 
pressure  value  and  the  procedure  is  repeated.  If  the  cntrojiics  do  not  agree, 
the  temperature  estimate  and  the  equilibrium  calculation  are  iterated  until 
the  entropies  do  agree. 


After  the  first  13  pressure  values  have  been  used  by  the  program, 
their  corresponding  equilibrium  and  thermal  values  for  the  STG  test  gas 
mixture  will  be  output.  Then  the  remaining  13  groups  of.  data  will  be 
calculated  and  printed.  Upon  completion  of  the  pressure  schedule,  the 
program  will  look  for  a  new  test  gas  mixture  in  the  input. 
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C  1.0000  0  1.0000  0.0  0.0  00  0.0  0.1400  M  0.0  G  298. ISO 
C  1  0000  0  2.0000  0.0  0.0  00  0.0  0.0250  M  O.O  G  298.150 
N  2*0000  0.0  0.0  0.0  00  0.0  0.2900  M  0.0  G  298.150 
AR  I.'OOOO  0.0  0.0  0.0  00  0.0  0.5450  M  0.0  G  298.150 
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